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ABSTRACT 

^^A  free-fall  oceanographic  instrument  has  been  used  to  measure 
vertical  microstructure  scale  gradients  of  horizontal  velocity,  tempera- 
ture and  electrical  conductivity.  The  velocity  gradients,  or  shears,  were 
measured  at  scales  between  3 and  40  cm  by  an  airfoil  shear  probe  whose 
specifications  and  calibration  procedure  are  discussed. 

Data  collected  in  the  equatorial  Atlantic  in  July  1974  indicated  a 
consistent  pattern  of  turbulence  near  the  velocity  core  of  the  Atlantic  Equa- 
torial Undercurrent.  (The  velocity  core  is  the  region  of  maximum  speed.  ) 
The  most  intense  turbulence  was  found  above  the  velocity  core  of  the  under- 
current. Turbulence  in  the  velocity  core  was  weak  and  intermittently 
spaced.  Below  the  core,  near  the  base  of  the  thermocline,  moderately  in- 
tense turbulence  was  found.  The  rate  of  viscous  dissipation  of  turbulent 

energy  has  been  estimated  from  the  shear  measurements,  and  typical 
-3  2 -3 

values  were  3x10  cm  sec  above  the  velocity  core. 

Spectra  of  the  shears  have  been  computed.  At  small  wavelengths 
the  measured  spectral  coefficients  fall  below  the  universal  Kolmogoroff 
spectrum.  This  discrepancy  between  the  two  spectra  is  attributed  to 
spatial  averaging  of  velocity  fluctuations  by  the  shear  probe.  The  estimates 
of  viscous  dissipation  include  a correction  for  this  spatial  averaging. 

An  energy  balance  has  been  determined  for  the  turbulent  velocity 
fluctviations.  Above  and  below  the  core  the  basic  balance  is  local  produc- 
tion of  turbulent  energy  equals  local  dissipation,  and  this  balance  gives  a 
v'ertical  eddy  viscosity  of  order  10  cm^sec”^ 


above  the  core.  The  equation 


iii 

of  the  energy  balance  of  the  average  motion  has  been  vertically  integrated 
at  the  equator  where  meridional  terms  are  assumed  small.  In  the  South 
Equatorial  Current  the  rate  of  energy  gain  from  the  average  zonal  wind 
stress  is  balanced  by  the  rate  of  energy  loss  to  the  zonal  pressure  gradient 
plus  the  rate  of  dissipation.  In  the  undercurrent,  above  the  core,  the  rate 
of  energy  gain  from  the  zonal  pressure  gradient  equals  the  rate  of  dissipa- 
tion within  the  uncertainty  of  the  measurements,  and  the  advection  term  is 
small  but  not  negligible. 
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INTRODUCTION 


The  energetics  of  equatorial  undercurrents  have  received  limited 
attention.  For  a flow  in  geostrophic  balance,  there  is  no  transfer  of 
energy  between  the  pressure  and  kinetic  energy  fields,  since  the  Coriolis 
force  can  do  no  work  and  flow  is  perpendicular  to  the  pressure  gradient. 

The  rate  of  energy  transfer  in  most  regions  of  the  oceans  depends  upon  the 
ageostrophic  flow  across  lines  of  constant  pressure,  which  is  relatively 
low  and  difficult  to  measure.  However,  at  the  equator  the  Coriolis  force 
disappears,  allowing  the  exchange  of  energy  between  the  pressure  and 
kinetic  energy  fields  to  be  more  easily  measured. 

Even  the  simplest  energy  balance  for  the  undercurrent  must  include 
the  work  done  by  the  zonal  pressure  gradient,  the  energy  advected  by  the 
mean  flow,  and  the  energy  dissipated  by  turbulence  (Wyrtki  and  Bennett, 
1’963).  The  rate  of  release  of  energy  from  the  zonal  pressure  gradient  can 
be  obtained  from  measurements  of  velocity  and  dynamic  height.  However, 
the  velocity  cores  of  the  equatorial  undercurrents  have  been  observed  to 
meander.  To  accurately  determine  the  velocity  field,  either  measurements 
at  several  locations  must  be  taken,  or  many  sections  across  the  equator 
must  be  sampled.  Until  1974,  such  observations  were  not  available.  The 
energy  advected  by  the  average  flow  remains  difficult  to  calculate  because 
of  limited  observations.  The  undercurrents  undergo  seasonal  changes,  as 
well  as  meanders  of  shorter  period.  Most  measurements  of  the  large 
scale  features  of  the  undercurrents  have  usually  been  at  a single  position 
and  extend  over  a relatively  short  time  interval,  so  it  has  been  difficult  to 


distinguish  temporal  variations  from  spatial  variations. 
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The  rate  of  dissipation  of  kinetic  energy  can  be  estimated  from 
velocity  measurements  at  scales  of  about  1 to  50  cm.  This  range  of  scales 
is  the  viscous  dissipation  range  at  which  kinetic  energy  is  converted  to 
heat  by  the  molecular  viscosity  of  the  fluid.  Instrumentation  to  sense  velo- 
city gradients  at  the  scales  and  magnitudes  found  in  the  undercurrents  has 
been  developed  only  recently.  Only  one  previous  attempt  has  been  made  to 
measure  these  gradients  in  the  undercurrent  (Belyaev  et  al,  1975b).  From 
measurements  of  temperature  fluctuations  and  the  use  of  theories  of  isotro- 
pic turbulence,  the  rate  of  dissipation  may  be  inferred,  as  was  done  by 
Williams  and  Gibson  (1974)  at  0®and  1®  N in  the  Pacific,  at  150®W,  at  a 
single  depth  at  each  latitude.  An  estimate  of  vertical  heat  transport  for 
certain  flows  can  be  made  from  measurements  of  the  mean  and  microscale 
temperature  gradients  (Osborn  and  Cox,  1972).  Such  measurements  in  the 
Pacific  at  the  equator  have  given  a qualitative  comparison  of  the  micro- 
scale velocity  gradients  and  the  intensity  of  turbulence  (Gregg,  1976). 
Although  theories  suggest  that  the  rate  of  dissipation  (£  ) should  be  a limit- 
ing factor  in  the  flow  of  the  undercurrent,  the  values  of  Z suggested  by 
these  studies  differ  from  the  value  derived  from  a study  by  Wyrtki  and 
Bennett  (1963)  based  upon  measurements  of  large  scale  features  in  the 
Pacific  Equatorial  Undercurrent. 

In  the  summer  of  1974,  participants  in  the  Global  Atmospheric 
Research  Program  undertook  a study  of  the  atmosphere  and  ocean  in  the 
tropical  Atlantic,  called  the  GARP  Atlantic  Tropical  Experiment  (GATE). 
During  this  project,  a study  of  the  microscale  velocity  gradients  in  and 


near  the  undercurrent  was  undertaken  with  instrumentation  developed  by 


T.R.  Osborn  and  myself.  The  study  was  conducted  from  the  R/V 
ATLANTIS  II  of  the  Woods  Hole  Oceanographic  Institute  during  June  and 
July  1974.  An  extensive  mapping  of  the  large  scale  velocities  and  dynamic 
heights  at  the  equator  was  undertaken  by  Dr.  Eli  Katz  and  Mr.  John  Bruce 
on  board  the  R/V  ATLANTIS  II,  and  by  numerous  investigators  on  board 
other  ships.  These  large  scale  measurements  make  it  possible  to  estimate 
some  of  the  transfers  of  the  average  kinetic  energy.  A comparison  between 
these  transfers  and  the  viscous  dissipation  of  turbulent  energy  is  possible 
for  a long  term  average  of  the  flow  at  the  equator. 

There  are  clearly  two  topics  here:  the  energetics  of  the  under- 
currents, and  of  microscale  oceanic  turbulence,  which  are  separated  not 
only  by  scale  size,  but  also  by  the  background  information  required  to 
understand  them.  In  the  next  section,  I will  discuss  these  topics  separate- 
ly, and  show  how  microscale  measurements  lead  to  a better  understanding 
of  the  overall  energy  balance  in  the  undercurrent. 

The  design  and  operation  of  the  instruments  will  be  discussed  in 
section  3,  as  well  as  the  procedure  in  the  field.  The  microstructure  of 
the  undercurrent  region,  the  nature  and  spectra  of  the  turbulence  there, 
and  the  energetics  of  the  undercurrent  will  be  described  in  section  4. 
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2.  BACKGROUND 

2.  1 Large  Scale  Features  of  the  Undercurrent 

There  have  been  numerous  models  constructed  to  explain  and  pre- 
dict the  nature  of  the  equatorial  undercurrents,  beginning  soon  after  the 
discovery  of  the  Pacific  Equatorial  Undercurrent,  as  described  by 
Cromwell,  Montgomery  and  Stroup  (1954).  A summary  of  such  models 
was  undertaken  by  Philander  (1973).  To  understand  why  equatorial 
currents  exist,  a simple  explanation  is  given  by  Neumann  and  Pierson 
(1966,  p.  453-454).  Figure  la  shows  a wind  and  current  pattern  typically 
found  near  the  equator.  The  doldrums  between  the  SE  and  NE  trade  winds 
are  located  near  B^north.  In  the  Southern  Hemisphere,  transport  in  the 
Ekman  layer  is  to  the  left  of  the  wind;  in  the  Northern  Hemisphere  this 
transport  is  to  the  right  of  the  wind,  causing  a convergence  of  the  surface 
waters  at  S^N  and  divergence  at  10®N  and  at  the  equator.  The  resulting 
-surface  slope  caused  by  such  convergences  and  divergences  is  illustrated 
in  figure  lb.  For  the  surface  currents  to  be  in  geostrophic  balance  with 
the  meridional  pressure  gradients  caused  by  these  surface  slopes,  they 


must  flow  westward  at  all  latitudes  except  between  5“and  10*N.  The  east- 
ward current  here  is  the  Equatorial  Countercurrent  in  the  Atlantic  Ocean. 
The  westward  current  present  to  the  south  of  S^N  is  the  South  Equatorial 
Current  and  that  north  of  1 O^N  is  the  North  Equatorial  Current.  This 
current  pattern  is  most  like  that  found  in  the  Atlantic  from  August  to 
November.  During  the  remainder  of  the  year,  when  the  winds  are  weaker, 
th('  countercurrent  is  frequently  absent.  In  the  Pacific,  this  countercurrent 


is  labellefl  the  North  Equatorial  Countercurrent  to  distinguish  it  from  the 


5 


North  Equatorial  Current 
Equatorial  Countercurrent 

South  Equatorial  Current 


N 


(b) 


'^iRuro  In  Possible  wind  distribution  in  equatorial  regions.  Long  ar- 
rows represent  wind  direction.  Short  arrows  sliow  possible 
F.knan  drift.  (from  Neumann  and  Pierson,  I'ihh  p.  452) 


Figiiro  lb  Idealized  topograpliy  of  the  sea  surface  across  the  equator 

caused  by  wind  pattern  of  figure  la.  Surface  water  divergen- 
ces or  convergences  are  denoted  by  PfV.  or  CONI',  and  direc- 
tion of  flow  is  <lenoted  by  L or  V.  (from  Neumann  and  Pier- 
son, Ihbb  p.  452) 


6 


South  Kquatoral  CkMinlu  riu  r runt  sonietinu-s  found  near  IO“soath. 

Whun  thu  winds  and  surfauu  l urrcnts  aru  present  as  shown  in  figu- 
res la  and  lb  the  undercurrent  is  set  up  in  the  followine  way:  the  west- 
ward surface  current  at  the  equator  causes  an  upward  surface  slope  to- 
ward the  west,  and  a downward  tilt  of  the  thermocline  toward  the  west. 
Associated  with  these  slopes  is  an  eastw-ard  pressure  force  in  and  above 
the  thermocline.  The  wind  stress  balances  or  even  exceeds  this  pressure 
force  in  the  surface  mixed  layer,  but  the  direct  wind  effects  cannot  pene- 
trate the  therniocUne.  Hence,  there  is  a westward  pressure  force  in  the 


thermocline  which  is  not  balanced  by  the  wind  stress.  Near  the  equator 

this  causes  an  eastward  flow  of  water  in  the  therniocline. 

"However,  wdth  increasing  distance  from  the  equator  the  effects 
of  the  Coriolis  force  become  more  and  more  significant  and 
tend  to  provide  a geostrophic  balance  of  the  zonal  pressure  gra- 
dient. This  seems  to  occur  very  close  to  the  equator,  probably 
as  close  as  for  1"  latitude.  Quantitative  evidence  of  approxima- 
te geostrophic  balance  within  the  undercurrent  to  latitudes  as 
low  as  i^has  been  presented  by  Montgomery  and  Stroup  (19b<H 
and  by  Metcalf  et  al  (1962).  Thus,  with  a downward  slope  of 
the  sc'a  surface  from  west  to  east  a little  north  of  the  equator, 
the  water  flows  with  a Southward  component,  and  a little 

south  of  the  equator  with  a Northward  component riiis 

leads  to  a water  transport  from  both  the  Northern  and  Southern 
Hemisphere  toward  the  equator  where  the  general  motion  is 
eastward  and  provides  ample  water  supply  lor  tlu'  F.ciuatorial 
Undercurrent."  (Neumann  and  Pierson  19f>6,  p.  ■4‘i4.) 

Cross-sections  of  temperature,  salinity  and  velocity  sampled  durinc 
the  ATLANTIS  IT  sections  along  24'’W,  ZH’W,  are  slun«.n  in  lieures  2 

to  4.  The  divergence  of  surface  currents  at  the  ec|uator.  cumliim-d  \\ith  in- 
flow from  the  sides  in  the  therniocline,  leads  to  vipwelliiu:  near  the  t-quator 
,ind  causes  the  ritlging  and  t roughing  of  the  isotherms  .shown  in  lienrc!,  2 to 


4.  This  structure  is  a charai  le  rist  ic  of  all  e(|ualorial  unde  tu  u r rent  s . 1 he 
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water  in  the  core  of  the  Atlantic  Equatorial  Undercurrent  derives  from  the 
North  Brazilian  Coastal  Current  which  has  high  salinity.  This  saline  core 
is  evident  in  these  cross-sections  and  has  been  found  in  this  undercurrent 
as  far  as  6“E,  where  the  undercurrent  was  observed  to  turn  away  from  the 
equator  (Rinkel,  et  al,  1966).  The  velocity  measurements  were  obtained 
from  an  over-the-side  current  meter,  and  are  relative  to  the  measurement 
at  300  m.  depth.  These  figures  are  taken  from  Bruce  and  Katz  (1976)  and 
are  based  on  data  they  collected  during  the  ATLANTIS  II  cruise. 


1.  2 Energy  Balance  for  the  Average  Flow 

It  is  a common  technique  for  incompressible  turbulent  flows,  first 
introduced  by  Osborne  Reynolds,  to  portion  the  velocity  and  pressure  fields 
into  average  and  fluctuating  components 
velocity  ij  ^ “^i  ’’’^i 

pressure  P = P + P*  2.  1 


where  P is  the  average  pressure  and  Uj^  is  the  average  velocity. 

When  these  values  are  included  in  the  energy  equations,  and  averaged,  one 
can  trace  the  energy  exchange  between  a pressure  field,  kinetic  energy  and 
heat.  In  theory,  one  uses  an  ensemble  average,  which  is  an  average  over 
all  realizations  of  the  flow.  In  practice  a time  or  space  average  is  measured 
within  the  flow  and  the  Ergodic  hypothesis  (Hinze,  1959,  p.  5)  is  invoked  so 
that  the  measurements  can  be  assumed  to  be  comparable  to  the  theory.  For 
time  averages  one  must  ask:  are  the  averages  over  a time  sufficiently  long 
to  include  all  or  most  of  the  fluctuations?  Some  periods  at  which  the 
Atlantic  Equatorial  Undercurrent  is  expected  to  contain  energy  arc: 


1 year  Annual  variations  in  the  winds  likely  produce  changes  in  the 

undercurrent  (Philander,  1973),  (Katz  et  al,  1976) 

16  days  Meandering  of  the  undercurrent  was  observed  in  the  summer 

f of  1974  during  GATE  (Duing  et  al,  1975) 

- 2 days 

3 riays  Pulsing  of  the  undercurrent  was  observed  during  GATE 

(Duing,  1975) 

12-25  Variations  due  to  tides,  with  the  semi-diurnal  tide  strongest 

hours  were  observed  during  GATE  (Rybnikov,  1975) 

several  min-  Internal  waves 

tes  to  several 

hours 

fractions  of  Turbulence 
seconds  to 
several  min- 
utes 

Prior  to  the  GATE  oceanographic  observations  in  1974,  there  were 
insufficient  observations  to  permit  distinction  between  variations  of  several 
days  and  annual  variations.  The  16  day  -Z  day  meander  appeared  in  the 
three  GATE  observ'ation  periods  of  21  days  each,  but  there  were  also  x’ari- 
ations  irom  one  observation  period  to  another  which  may  or  may  not  be 
part  of  an  annual  variation.  The  observations  during  the  GATE  project 
form  the  longest  set  available  and  in  many  cases  permit  averaging  over  the 
16  day  wave.  No  averages  over  longer  variations  are  available.  The  effect 
of  short  sampling  times  upon  the  observations  of  the  mean  gradients  may  be 
large  in  some  cases. 

To  obtain  the  energy  equation  of  the  average  flow,  the  average  and 
fluctuating  values  of  velocity  and  pressure  are  substituted  into  the  Navier- 
Stokes  equation  for  a fluid  of  constant  density  and  \ iscosity.  The  coordina- 
te system  is  chosen  with  Xj  eastward  along  the  eejuator,  X7  northwa 
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oriuin  at  the  equator,  and  positive  upward  and  zero  at  a level  surface 
chosen  to  be  close  to  the  actual  ocean  surface.  Each  term  in  the  Navier- 
Stokes  equation  is  multiplied  by  Ui  and  an  ensemble  average  of  each  term 
is  taken.  When  all  the  terms  are  summed,  the  rate  of  change  of  kinetic 
energy  per  unit  mass  and  per  unit  time  is,  in  tensor  form  (liinze,  1959,  p. 

Q^)  -c,U^ 

^xj  f dxjalij  ^ -Xj  ^ 

where  the  ox-erbar  denotes  an  ensemble  average. 


=_L[UNU2-U3l 

z 2 L J 


2.  3 


is  the  kinetic  energy  per  unit  mass,  and  use  has  been  made  of  the  continu- 
ity equation  for  incompressible  flows 


= o 

a:; 


2.  4 


The  term  -gU^  has  been  added  to  Hinze's  equations  because,  while  he  is 
able  to  use  a reduced  pressure  (total  minus  hydrostatic),  such  an  approach 
cannot  be  used  in  this  case  because  of  the  presence  of  a free  sloping  sur- 
face'. 

It  is  assumed  that  the  hydrostatic  balance  is  valid,  and  hence 


Us  ^ P ^ "<9  U5 
f ^ 


2.  5 


To  reduce  the  number  of  terms,  the  -gU3  term  may  be  omitted  by  adopting 

the  convention  that  the  U-^  term  is  to  be  summed  over  i=l  and  2 but  not  3. 

^ 3xi 

There  is  no  Coriolis  term  in  equation  2.2.  Although  the  Coriolis 


fon  e can  alte'r  the  direction  of  flow,  it  cannot  change  the  speed  of  flow 


1 3 


and  cannot  alter  the  magnitude  of  the  kinetic  energy. 

If  a small  density  variation  is  present,  the  density  field  can  be 

portioned  into  an  average  and  a fluctuating  component  f ^ + p'  in  the 

Navier-Stokes  equations.  K ^ , the  Bous sinesq  approximation  may  be 

invoked  to  simplify  the  equations.  As  outlined  by  van  Mieghem  (1973),  if 

^ < 10”  , the  Boussinesq  approximation  holds:  variations  in  the  density 

of  the  fluid  are  only  important  where  they  affect  the  buoyancy.  For  the 

undercurrent,  < 2x10”  . Hence,  in  equation  2.2,  the  density  p may  be 

? “ 

replaced  by  the  average  density^,  but  it  must  be  outside  the  derivative  for 
the  equation  to  be  correct. 

The  term  in  equation  2.  2 is  the  sum  of  the  rate  of  dissipa- 

tion  per  unit  mass  of  the  kinetic  energy,  and  the  rate  of  work  done  per  unit 
mass  by  the  viscous  stresses.  For  oceanic  flows,  this  term  is  very  small 
compared  to  others  and  may  be  neglected. 

The  last  term  in  equation  2.  2 may  be  rewritten 


- Ui  'iUj 


“ ^ LX'. 

Xj  ^ Xj 


2.  6 


Equation  2.  2 is  now 

~ Ui (y  - ^ U ; Hi  U j -f  u.UjC^U; 

f 5xi  ^ -Xj 

I II  Ilia  mb  IV 

For  the  meanings  of  the  individual  terms  in  equation  2.  7 wc  have: 

I The  local  rate  of  change  of  average  kinetic  energy  per  vinit  mass, 

for  example  a seasonal  change. 


II 


The  rate  of  work  done  per  unit  mass  by  the  pressure  field.  Mea- 
surements of  the  anomaly  of  dynamic  height  and  the  average  velo- 
city along  the  equator  permit  the  evaluation  of  this  term. 
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- 

Ilia 

Illb 

Divergence  of  transport  of  kinetic  energy  per  unit  mass  by  the 
average  flow  and  the  turbulent  diffusion  of  average  energy. 

- 

IV 

The  rate  of  exchange  of  energy  between  the  kinetic  energy  of  the 
fluctuating  motion  and  the  kinetic  energy  of  the  average  motion. 
It  is  the  product  of  the  Reynolds  stress  per  unit  mass  u,  Uj 
and  the  average  shear, 

A distinction  has  been  made  by  Bowden  (1962)  between  fluctuations 


due  to  turbulence  and  those  due  to  two  dimensional  eddies.  Turbulence  is 
a term  for  random  fluctuations  of  all  three  components  of  velocity,  requir- 
ing a transfer  of  energy  from  large  scales  of  motion  to  smaller  scales  and 
finally,  at  the  smallest  scales,  the  transfer  of  energy  to  heat.  Tvirbulence 
is  always  a sink  for  kinetic  energy  of  the  mean  flow,  and  causes  mixing. 

Large  scale  two  dimensional  eddies  (or  meanders)  are  often  found  in  the 
ocean  or  atmosphere,  may  eitlier  extract  energy  from  or  contribute  enerpy 
to  the  mean  flow  and  need  not  cause  vertical  mixing.  The  contribution  of 
both  of  these  interactions  to  the  kinetic  energy  of  the  average  flow  is  repre- 

1 
I 

sented  by  term  IV. 

Estimates  of  some  of  the  terms  in  equation  2.  7 are  available  from 
oijservations  during  GATE,  and  will  be  discussed  in  section  4.  Term  I may 
be  estimated,  as  well  as  the  zonal  components  of  terms  II  and  Ilia.  Regions 
can  be  chosen  over  which  the  contributions  of  other  divergence  terms  are 
small.  How’ever,  there  are  no  measurements  in  equatorial  regions  of 
term  IV  in  equation  2.7.  It  is  possible,  however,  to  evaluate  some  terms 
in  the  kinetic  energy  equation  for  the  fluctuating  flow.  The  term  ic. u ; 

^ ^,1 

appears  in  this  equation  as  a sovirce  of  kinetic  energy  of  turbulent  fluctua- 
tions, anrl  as  a source  or  sink  of  kinetic  energy  of  the  two-dimensional 


eflflies.  If  other  terms  in  the  efjuation  are  measured  or  small,  the  term 
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may  be  deduced. 

2.  3 Energy  Balance  for  the  Fluctuating  Flow 

The  equation  for  the  kinetic  energy  of  a constant  density,  fluctuating 
flow  can  be  derived  in  a way  similar  to  that  of  the  equation  for  the  average 
flow.  The  Navier-Stokes  equation  is  multiplied  by  U.,  averaged  and  summed. 
The  equation  in  tensor  form  is  (Hinze,  1959,  p.  65): 

JT  5^  f cJ^C;  'Sxj  2.8 

I II  III  IV  V 

where 

2.9 

is  the  kinetic  energy  per  unit  mass  of  the  fluctuating  motion. 

Equation  2.  8 states  that  the  rate  of  change  (I)  of  kinetic  energy  of 

the  fluctuating  motion  per  unit  mass  of  the  fluid  following  the  mean  motion 

(i.c.  , ) is  equal  to  (II)  the  convective  diffusion  by  the  fluctua- 
nt ^Xj 

ting  motion  of  the  total  energy  of  the  fluctuating  motion,  plus  (III)  the  energy 
transferred  to  or  from  the  mean  motion  through  the  stresses  of  the  fluctuat- 
ing motion,  plus  (IV)  the  work  done  per  unit  mass  and  time  by  the  viscous 
stresses  of  the  fluctuating  motion,  plus  (V)  the  dissipation  per  unit  mass  of 
th('  fluctuating  motion  (Hinze,  1959).  As  shown  by  Tennekcs  and  Lumloy 
(1972,  p.  70),  term  IV  is  of  order  (Reynolds  number)*  * times  term  V and 
may  be  neglected. 

The  last  term  is  frequently  written 


f 

\ 

= -J2l 

< ^ U : + \ 

= - € 

2.  10 

^x, 

1 

2 ' 

^dXj 

16 


which  shows  that  the  viscous  dissipation  term  in  equation  2.  8 is  always 
negative  definite,  and  alw'ays  decreases  the  value  of  converting  the  kine- 
tic energy  into  heat. 

The  presence  of  small  density  variations  alters  the  form  of  equation 
2.  8.  As  was  noted  for  the  equation  of  the  kinetic  energy  of  the  average  flow, 
the  Houssinesq  approximation  is  valid  for  the  undercurrent.  When  the  den- 
sity field  is  portioned  into  mean  and  fluctuating  components,  and  the 
Boussinesq  approximation  applied,  the  density  p may  be  replaced  by  the 
a\crage  density  p in  equation  2.8,  (however,  it  must  be  outside  the  deriva- 
tive) and  an  extra  term  -p'uj(^/p  is  added  to  the  right  side  of  the  equation 
(van  Mieghem,  1973). 

Ir>  the  presence  of  density  variations,  is  the  buoyancy  force 

per  unit  mass  acting  on  a fluid  element,  and  - is  the  average  work  done 
per  unit  mass  by  the  buoyancy  force.  In  a stably  stratified  fluid,  this  term 
is  negative,  indicating  that  kinetic  energy  is  lost  to  the  potential  energy 
field.  When  the  stratification  is  unstable,  kinetic  energy  may  be  gained 
from  the  potential  energy  field,  and  this  term  is  positive  (Monin  and  Yaglom, 
1971,  p.  397). 

The  equation  for  the  kinetic  energy  per  unit  mass  of  the  fluctuating 
velocity  is  now 

= -I  - u,Uj(^Ui  -£ 

dt  P ^ 2.  1 1 

I II  III  IV  V VI 

The  fluctuating  components  considered  hero  include  the  meanders, 
tides,  internal  waves,  and  turbulence.  The  expected  effect  of  the  meanders 
lias  been  fliscussed;  the  tides  and  internal  waves  can  transmit  energy  from 


one  region  to  another,  and  can  carry  energy  into  and  out  of  a region.  How- 


ever, unless  they  interact  with  turbulence  or  generate  instabilities  they  are 
not  very  dissipative,  and  do  little  mixing  (for  a wave  an  excellent  approxima- 
tion is  Ll:uj  = 0 for  i ^ j);  hence  their  contributions  to  terms  V and  VI  are 
small.  To  produce  large  values  of  e , turbulence  is  required.  The  evalua- 
tion of  terms  II  and  III  is  difficult,  and  has  not  been  attempted  for  oceanic 
turbulence.  There  are  no  measurements  of  term  IV,  as  mentioned  in  the 
previous  section,  and  no  measurements  for  term  VI. 

With  such  a dearth  of  measurements,  one  is  forced  to  some  assump- 
tions. In  the  Introduction  it  was  noted  that  the  viscous  dissipation  of  turbu- 
lent energy  was  an  important  sink  for  the  kinetic  energy  of  the  average  flow. 
This  energy  is  passed  from  the  average  flow  to  the  turbulence  through  term 
IV,  which  appears  in  equation  2.  7 and  equation  2.  1 1 . There  are  no  equator- 
ial measurements  of  term  IV,  but  I will  show  in  section  4 that  it  is  similar 
to  e in  regions  of  relatively  intense  turbulence  at  the  equator,  and  e can  be 
estimated  from  the  velocity  microstructure  measurements. 

To  determine  the  dissipation,  the  degree  of  isotropy  must  be  known 
because  it  is  not  possible  to  measure  all  the  terms  in  the  expression  for  c 
(equation  2.  10).  In  flows  of  large  Reynolds  numbers,  there  may  be  several 
orders  of  magnitude  difference  in  the  size  of  scales  at  which  turbulent 
energy  is  generated  from  the  mean  flow  and  the  scales  at  which  it  is  conver- 
ted to  heat  by  the  viscous  forces.  In  such  cases  Kolmogoroff  first  proposed 
in  1941  that  there  is  a tendency  for  flows  at  small  scales  where  viscous 
dissipation  occurs  to  become  isotropic,  a condition  termed  local  isotropy. 


In  this  case  (Monin  and  Yaglom,  1975,  p.  353-354): 
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O ' ' 


1.  12a 


Thus  measurements  of  a single  gradient  at  small  scales  will  give  the  full 
value  of  the  dissipation.  Writing  out  €.  in  full  from  equation  2,  10  we  have 


j_  ^ ^ Li j 

2 ^ ^ 

lax?/ 

+ au;  f'  +fSU2.  + ^U3 

la^i  ax,/  vjxj  axW  ax?/ 


2.  12b 


In  the  isotropic  case 


/ <5Ui 

= V 

^ ^x,  / 

( ^xz/ 

15^/ 

and  terms  like 


/,^u,  -f 

v ^Xj/ 

^ ax,  / 

and  cross  terms  like 

Z =z  — ( ^ 

dx,  axj  ” ^ J 

The  instrumentation  we  used  at  the  equator  (and  described  in  section  3) 
measures  ith  and  directly  (although  some  corrections  must  be  made 

of  small  scale  resolution  limits  of  the  probe).  Thus  in  the  isotropic  case, 
the  sum  of  the  quantities  measured  isj|£  . Isotropy  is  needed  only  when  we 
convert  what  we  measure  to  an  estimate  of  the  total  dissipation.  Note  that 
even  if  the  turbulence  is  anisotropic  we  are  measuring  directly  terms  which 
make  a large  contribution  tot  , since  equation  2.  10  is  a general  relation 
not  di'pcndent  on  isotropy. 
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2.  4 Spectral  Representation  of  Turbulence 

Some  checks  of  local  isotropy  may  be  made  using  the  spectral  rep- 
resentation of  turbulence.  G.  I.  Taylor  first  introduced  the  energy  density 
function  to  describe  the  concept  of  energy  in  turbulence  associated  with  a 
scale  size.  I will  use  the  notation  E(k)  to  denote  the  three  dimensional 
energy  density  function;  E(k)dk  denotes  the  energy  of  all  fluctuations  having 
a w'avenumber  magnitude  between  k and  k+dk.  The  one  dimensional  spectr- 
al function  cp^  ^ (kj)  denotes  the  energy  density  of  the  Uj  component  of  velo- 
city associated  with  all  fluctuations  having  a wavenumber  component  k^  in 
the  Xj  direction,  and(P22  the  energy  density  of  the  Up  component  of 

velocity  associated  with  all  fluctuations  having  a wavenumber  component 
kj  in  the  x j direction.  The  three  velocities  are  related  to  the  spectral 
functions  by 
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In  a locally  isotropic  region  of  wavenumber  space  one  finds  (Monin  and 
Yaglom,  1975,  p.  54): 
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and  also 

- O 2.19 

where  ) is  the  spectrum  of  JTu^  associated  with  a wavenumber  kj  such 

oO 

f (p,^ik,)dk,  = uuj 


that 
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Equations  2.  17  to  2.  19  are  necessary  (but  not  sufficient)  conditions  for  the 
presence  of  isotropy  at  the  scale  size  corresponding  to  the  wavenumber  kj. 
Theories  of  isotropic  turbulence  are  generally  expressed  in  terms  of  E(k) 
which  cannot  be  measured.  Measurements  of  turbulence  usually  sense  com- 
ponents leading  to  estimates  of  j (k^),  Cp  22  ‘Pl2^^1^’  which  are 

compared  with  theoretical  relationships  through  equations  2.  17  to  2.  19. 

The  viscous  dissipation  for  locally  isotropic  turbulence  is 

C - Zv  ^ £.(k)dk 

O 

or,  in  terms  of  measured  values  of  ^ (k^)  and  is  (Monin  and 

Yaclom,  1975,  p.  405): 


= (p^/kJdk. 
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The  spectral  representation  can  be  used  to  check  for  local  isotropy 
within  the  dissipation  region  and  to  determine  the  dissipation.  A second 
arlvantage  is  that  if  the  turbulence  is  locally  isotropic,  the  form  of  the 
spectral  density  functions  can  be  predicted,  and  conversely,  if  the  shape 
of  one  of  the  Cp,|  (k  j ),  <p2Z^'^P’  <p33^’^l  ' measured,  then  assuming 

isotropy,  € can  be  rlet<'rmined  without  measurements  at  dissipation  scales. 
These  determinations  rely  upon  the  two  Kolmogoroff  hypotheses. 
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For  large  Reynolds  numbers,  Kolmogoroff  proposed  that  the  statis- 
tical properties  of  the  flow  at  small  scales  are  independent  of  the  large 
scale  flow,  except  for  the  average  rate  e at  which  energy  is  passed  from 
the  large  scales  to  the  smaller.  These  statistical  properties  will  also  de- 
tjend  upon  the  kinematic  viscosity  v>  and  the  wavenumber  k,  A characteris- 
tic wavenumber  can  be  formed  from  the  viscosity  and  the  dissipation 

ks  - Cl  2..  23 


where  k^  is  the  Kolmogoroff  wavenumber.  Dimensional  analysis  gives. 


(p.ck.) 
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Both  F(^/k_)  and  F-,(^/k-)  are  functions  of  the  universal  non-dimensional 
wavenumber  (k/kg).  (The  subscript  is  traditionally  omitted  from  the  wave 
number  kj  within  the  brackets  but  it  is  understood  that  this  is  kj  , not  the 
magnitude  of  the  vector  wavenumber.  ) The  dissipation  spectra  are 


where 

G ( %5 ) ^ ) 


and 


are  also  universal  non-dimensional  spectra. 
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For  yet  larger  Reynolds  numbers,  the  second  Kolmogoroff  hypothe- 
sis proposes  that  only  e and  k are  the  relevant  parameters  within  a locally 
isotropic  inertial  subrange.  Scales  within  this  subrange  are  too  small  to 
be  directly  influenced  by  the  mean  flow,  and  too  large  to  be  influenced  by 
visi  osity.  Dimensional  analysis  gives 

(p„(k,)  = k £ ’ fe,  and  using  2.  18  2.  29 

3 ''  2.30 

where  K'  is  a universal  constant.  A summary  by  Monin  and  Yaglom  (1975) 
of  experiments  gives  = 0.  5 as  the  most  reasonable  value  suggested  by 
experimental  data.  Oceanic  work  by  Nasmyth  (personal  communication) 
gives  K'  = 0.52.  The  log  F(^/kg)  curve  derived  by  Nasmyth  (personal 
communication)  from  a re-evaluation  of  his  work  reported  earlier  (Nasmyth, 
1970)  is  shown  in  figure  5.  The  log  F2(  /kg)  curve,  calculated  from  data 
for  the  log  F(^/ks)  curve  using  equation  2.  24  is  shown  alongside.  Values 
of  log  Ci(k/kg)  and  G2(k/kg)  h ave  been  calculated  from  the  log  F(k/kg) 
i.  urve  using  equations  2.  27  and  2.  28  and  are  also  plotted  in  figure  5. 

To  determines  accurately  from  spectra  the  turbulence  must  be  local- 
ly isotropic  at  least  through  the  dissipation  portion  of  the  spectrum,  and  the 
Kolmogoroff  hypothesc-s  niust  be  followed  for  the  portion  of  the  spectrum 
scnserl.  It  is  known  that  the  Cp^  j(k^)  spectrum  often  follow'S  a -5/3  power 
dependt'nce  upon  k^  at  wavenumbers  where  the  turbulence  is  anisotropic 
(e.g.  , Wi'iler  anfl  Burling,  19f>7)  and  care  must  be  taken  in  analyzing  the 


ri-.-vilts.  It  the  fluid  is  stratifii'd,  buoyancy  effects  may  influenci'  thi-  spectra. 
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5 Buoyancy  Effects  on  Turbulence 

The  discussion  of  the  spectral  representation  of  turbulence  until 
now  has  assumed  that  effects  of  buoyancy  are  small.  When  a fluid  is  stably 
stratified,  the  turbulence  can  be  altered  in  several  ways.  To  extract  tur- 
bulent energy  from  the  average  flow,  the  Reynolds  stress  -^u^  u^  must  be 

appreciable.  However,  turbulent  energy  permit  mass  contained  within  the 

2 

vertical  velocity,  (u,  is  not  prodviced  directly  from  the  average  horizon- 
"2 

tal  flow,  but  must  be  transferred  from  other  turbulence  components  by  the 
]3ressure  forces.  The  buoyancy  sink  -u^p'g/p  extracts  energy  directly  from 
the  vertical  component.  Stewart  (1959)  suggests  that  the  extraction  processes 
are  more  efficient  than  the  production  processes  for  the  vertical  component 
of  iurlnilence,  because  the  production  is  indirect  and  must  involve  viscous 
losses  in  other  components  along  the  way.  To  maintain  turbulence,  the 
fhix  Richardson  number  (the  ratio  of  buoyancy  sink  to  total  shear  production) 
must  be  significantly  less  than  one  (Pond,  1973).  Turner  (1973)  indicates 
tliat  a flux  Richardson  number  greater  than  0.  15  may  be*  sufficient  to  elimin- 
ate turbulence.  Because  the  buoyancy  sink  removes  energy  only  from  the 
u^  component,  the  turbulence  may  not  be  locally  isotropic  over  the  same 
ranee  of  wavenumbers  as  for  the  non-stratified  flow,  l.ocal  isotropy  may 
only  be  found  at  higher  wavenumbers. 

Internal  waves  are  possible  with  frcciuencies  up  to  the  Brunt-\ais- 


ala  frefiuency. 
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These  waves  may  contribute  energy  to  the  . ‘P  >2  ^”'*^33 
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low  frequencies.  These  effects  are  discussed  by  Stewart  (1969).  Internal 
waves  in  a laminar  flow  can  interact  to  generate  patches  of  turbulence. 
Patches  of  turbulence  can  lose  energy  to  internal  waves  which  radiate 
energy  away. 

There  is  some  evidence  that  for  turbulence  of  sufficiently  high 
Reynolds  numbers,  stable  stratification  influences  the  energy  density 
spectra  only  at  wavelengths  greater  than  some  value  lb.  Turbulent  fluctua- 
tions of  wavelengths  less  than  iy  may  not  feel  the  large  scale  stratification, 
and  may  be  locally  isotropic.  In  a stably  stratified  atmospheric  boundary 
layer,  experiments  have  found  evidence  consistent  with  local  isotropy  for 
\'a lues  of  f L/ll ~ 1 0 ( ^ = 0.  1 L)  (Kaimal  et  al,  1972)  where  L is  the  Monin- 
Obukhov  length  and  A is  the  wavelength.  Away  from  the  surface  in  the 
atmosphere'  or  ocean,  the  length  L is  no  longer  the  relevant  length  scale, 
and  various  other  scales  have  been  proposed.  Ozmidov  (1965)  has  proposed 
the  criterion 

A = 

2.  32a 

and  Obukhov  (1959)  has  suggested  the  length  scale  for  regions  where  temp- 
(-rature  fluctuations  dominate  tlie  density  fluctuations  should  bt' 

where  (S> is  tlu'  rate  of  dissipation  of  temperature  f hictuat  ions , k-t 
is  the  moh’cular  diffusivity  for  heat  and  T^  is  the  average  temperature  in 
flc'grees  Kelvin.  Utifortunately  there  have  been  no  rigorous  tests  of  iso- 
tropy using  the  criteria  of  ef|uations  2.  17  to  2.  19,  although  several  inves- 

-5/3 


tigators  have  obsc-rvefl  a k 


power  law  dependence'  for  Cp  | | ( 


J 
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iP^^(k^):  (Shur,  1962,  Vinnichenko,  1969,  for  clear  air  turbulence  in  the 

atmosphere;  Seitz,  1973,  for  turbulence  in  an  estuary;  Belyaev,  et  al, 

1975a,  for  turbulence  in  the  ocean).  Measurements  by  Nasmyth  (1970)  of 

patches  of  turbulence  in  the  winter  thermocline  in  the  North  Pacific  Ocean 

revealed  that  the  spectra  did  not  conform  to  a universal  Kolmogoroff 

spectr\ini  even  for  scales  in  the  viscous  dissipation  region. 

If  the  turbulence  in  the  dissipation  range  is  anisotropic  due  to  stable 

density  gradients,  u^  may  not  be  suppressed  completely,  or  there  would  be 

no  turbulence.  Idowever.  if  u^  tends  to  zero,  but  the  isotropic  relations  hold 

for  u , and  u_,  then  £=1:  M ; that  is,  the  true  dissipation  would  be  ^ ^ / 1 5 

12  ^*-1' 

of  the  value  estimated  from  the  measurements  using  the  isotropic  relation. 
The  situation  cannot  be  this  extreme,  although  in  the  anisotropic  case 

<y  X j 

and  may  become  somewhat  larger  relative  to  terms  like 

etc.  It  appears  that  an  estimate  of  the  rate  of  dissipation  based  upon  our 
measurements  of  and  and  assuming  isotropy  may  overestimate  c 

ot)C.3 

by  almost  50"o  in  the  most  extreme  case,  but  likely  the  error  is  smaller. 


27 


2.  6 Observations  of  Oceanic  Turbulence 

To  measure  any  of  the  terms  in  the  expression  for  £ (equation  2.  10) 

velocity  gradients  at  scale  sizes  at  which  viscous  dissipation  takes  place 

-3  - 

must  be  resolved.  The  Kolmogoroff  wavenumber  = (£  V and  the 

corresponding  Kolmogoroff  wavelength  A g =^'''^/kg  characterize  the  wave- 

nunibers  and  wavelengths  at  which  the  rate  of  viscous  dissipation  is  high.  I 

For  turbulence  which  is  locally  isotropic  at  dissipation  scales,  the  peak  of  ; 

the  dissipation  spectrum  G(k/kg)  is  near  a wavenumber  k=0.  1kg.  (Pond,  s 

1965.  ) Oceanic  turbulence  has  a dissipation  spectrum  extending  from 
small  scales  of  1 to  6 cm  (Gregg,  1973),  up  to  about  1 metre,  and  velo- 
city fluctuations  at  these  wavelengths  are  referred  to  as  velocity  micro- 
structure. The  intensity  of  the  microscale  velocity  gradients  can  often  be 

used  as  a q 'alitative  indicator  of  the  intensity  of  turbulence  (i.e.  , the  inten-  \ 

sity  of  the  large  scale  energy  containing  fluctuations)  because  the  small 
scale  fluctuations  can  adapt  quickly  to  changes  of  the  large  scales,  and 

increase  in  amplitude  as  the  large  scale  fluctuations  increase  in  amplitude.  1 

Observations  by  Grant  et  al.(1962)  in  a tidal  channel  in  British 
Columbia  revealed  very  active  turbulence  by  using  a hot  film  probe  towed  ; 

behind  a ship.  They  measured  values  of  € from  0.  01  to  1 cm'^sec  ' and 
found  an  extensiv'e  minus  five  thirds  region  in  the  spectrum  of  (p^^lk^). 

Reynolds  numbers  in  the  channel  are  as  high  as  3x10  . I.ater  observations 
by  Grant  et  al.  (196k)  in  the  open  ocean  show  that  in  the  upper  mi.xed  layer 
measurable  turb<ilence  is  almost  always  foiind.  In  and  below  (he  thermo- 

( line,  thj  fraction  of  the  water  which  has  nu'a-surable  turbulence  is  small,  I 

oftf’n  less  than  10%.  Average  values  of  £ in  the  upper  mixed  layer  were  . 
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0.  025  cni‘"sec  , and  about  two  orders  of  magnitude  less  than  90  metres 
depth. 

Nasmyth  (1970)  obtained  measurements  of  turbulence  in  the  winter 

thermocline  in  the  North  Pacific  which  reveal  patchiness  of  the  turbulence. 

hU'  towed  the  hot  film  probe  used  by  Grant  et  al,  (1962,  1968)  through  the 

ocean,  but  added  a cycling  mode,  whereby  the  probe  could  be  alternately 

raised  and  lowered  to  sample  along  a sawtooth  profile.  Values  of  C ranged 
- 2 -4  2 - 3 

from  about  2x10  to  2x10  cm  sec  within  the  patches  of  turbulence, 
v'.hieh  occupied  about  15%  of  the  observed  path  during  one  set  of  observations, 
and  3 0%  during  another. 

A free  fall  instrument  used  by  Osborn  (1974)  to  sense  the  (k^) 

and  (k^)  components  of  velocity  microstructurc  in  Fiowe  Sound,  British 

Columbia,  revealed  large  variations  in  €.  over  short  distances.  An  examina- 
tion of  a typical  5 metre  thick  patch  showed  estimated  dissipation  there  of 

-4  2 - 3 

about  3x10  cm  sec 

Measurements  of  oceanic  turbulence  by  Belyaev  et  al  (1975  a,  b) 
have  revealed  large  changes  in  turbulent  intensity  in  the  horizontal  and  ver- 
tit  al,  with  the  magnitude  of  i.  varying  by  a factor  of  100. 

Observations  by  Woods  and  Fosberry  (1969)  in  the  seasonal  thermo- 
cline in  the  Mediterranean  have  revealed  that  shear  instabilities  can  grow 
on  a troiieh  f)r  crc'St  of  an  internal  wave  ])ropagating  through  a region. 

Slu'ar  instabilities  grow  from  small  disturbances  which  are  able  to  extract 
|■nl■|■L'y  from  a larger  scab'  shear,  roll  up  into  billows,  i-ventually  lumome 
g r,iv  itat  ional  ly  unstable  and  giuierate  turbulencu'.  'I'hese  i nvi-stigat  ors  noti'd 
that  the  instabilities  grew  on  regions  uhere  tlu-  vertiial  density  grarlient 


was  large,  called  "sheets".  The  shear  instabilities  generated  turbulence 
in  patches,  mixed  the  dye  over  a depth  of  up  to  20  cm,  and  were  smeared 
out  along  the  sheet  by  the  large  scale  shear. 

Theoretical  studies  by  Miles  and  Howard  (1964)  have  shown  that 
small  distvirbances  will  not  grow  in  a laminar  flow  if  the  gradient  Richard- 
son number 


Ri 


-g  ^ 
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is  greater  than  The  observations  of  Woods  and  Fosberry  (1969)  suggest 
that  in  and  below  the  thermocline  the  average  shear  is  too  small  to  main- 
tain turbulence,  but  that  internal  waves  passing  through  can  increase  the 
shear  in  a region  long  enough  to  allow  shear  instabilities  to  grow  and  be- 
come turbulent.  Measurements  in  the  Atlantic  by  Sanford  (1975)  have 
shown  that  10  m average  values  of  Ri^  below  the  thermocline  are  greater 
than,  but  often  close  to  0.25.  However,  Sanford  and  also  .Stewart  (1969) 
suggest  that  if  there  are  non-uniform  density  and  velocity  gradients  within 
the  region,  then  the  gradient  Richardson  number  must  be  smaller  for  a 
portion  of  the  region.  For  example,  if  within  a 10  metre  layer,  the  velo- 
city and  density  gradients  were  concentrated  within  a few  metre-s,  the 
actual  Richardson  numl)er  would  be  much  lower  within  these'  few  nic'tri'S 
than  indicated  by  a 10  metre  average  gradient. 

Thorpe  (1973)  found  shear  instabil  it  ic'S  would  grow  and  develop  into 
turl)ulence  in  a tank  only  for  \alues  of  Ri^,  much  h'ss  than  0.25.  Uecc'nt 
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laboratory  studies  by  Koop  (1976)  have  shown  shear  instabilities  to  grow  into 

billows  for  values  of  Ri  less  than  about  0.  07.  For  larger  values,  waves 

g 

grew  and  their  crests  broke  off  and  mixed  into  the  flow.  Orlanski  and 
Bryan  (1969)  have  proposed  that  internal  waves  may  grow  unstable  and  over- 
turn much  like  plunging  breakers  on  a beach.  Garrett  and  Munk  (1972)  con- 
clude that  shear  instabilities  growing  on  internal  waves  area  more  frequent 
source  of  turbulence  than  overturning  internal  waves,  but  they  point  out  that 
a single  overturning  internal  wave  may  generate  much  more  turbulent  kine- 
tic energy  than  a single  shear  instability. 

Double  diffusion  instabilities  extract  energy  from  the  potential  ener- 
gy of  the  fluid  to  drive  convective  regimes.  These  have  been  summarized 
by  Turner  (1973)  and  Gregg  (1973,  1975),  When  warm  salty  water  overlies 
colder  fresher  water,  salt  fingering  may  be  found.  Williams  (1974)  obser- 
ved such  fingering  below  the  Mediterranean  outflow  in  the  Atlantic,  of 
scales  about  6 mm  in  diameter  and  24  cm  high.  When  cold  fresh  water  over- 
lies  warm  salty  water,  layering  may  be  found.  Such  layers  have  been  obser- 
ved for  example  in  the  Arctic  by  Neshyba  and  Neal  (1971)  and  in  several 
other  oceanic  regions.  Thicknesses  of  the  layers  vary  from  less  than  a 
metre  to  up  to  20  metres  (Gregg,  1975).  Although  double  diffusion  regimes 
have  been  positively  identified  only  in  few  regions,  they  may  influence  tur- 
bulence in  much  of  the  ocean.  Gargett  (1976)  has  noted  that  there  is  fairly 
strong  statistical  evidence  for  increased  microstructure  activity  in  regions 
wlu-re  the  local  vertical  gradients  of  temperature  and  salinity  are  suitable 
for  flouble  fliffusion  processes. 
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?,.  7 Observation  in  the  Undercurrents 

Measurements  by  Jones  (1973)  in  the  Pacific  Equatorial  Undercurrent 
show  low  values  of  the  difference  Richardson  number  (a  finite  difference 
approximation  for  Rig) 

ft,,  = -1  2.34 

Values  of  Ri^  were  less  than  one  quarter  over  5 to  1 0 metre  depth  intervals 
between  196  and  272  metres  within  the  thermostad,  which  is  a layer  having 
a temperature  from  11.5  to  13.  5®C,  located  from  about  100  to  300  metres 
below  the  surface.  This  layer  is  characterized  by  small  changes  in  density, 
temperature  and  salinity  with  depth.  It  appears  to  increase  in  thickness, 
temperature  and  salinity  as  it  flows  from  west  to  east,  and  may  be  formed 
by  vertical  or  horizontal  mixing.  Jones  feels  that  the  low  values  of  Ri^ 
Indicate  that  vertical  mixing  should  be  more  important. 

Williams  and  Gibson  (1974)  have  measured  temperature  fluctiiations 
at  microstructure  scales  in  the  Pacific  Equatorial  Undercurrent.  From  the 
wavenumber  at  which  the  <p.prj,(kj ) vs  log  kj  spectrum  acquired  a -1 

2 -3 

.slope,  they  calculate  a viscous  dissipation  of  0.08  cm  sec  at  the  Equator 
and  at  1®  North.  These  calculations  require  that  conditions  for  Batchelor's 
(1959)  theory  for  spectra  of  temperature  fluctuations  be  met.  However, 
studies  by  Gregg  (1976)  in  the  Pacific  Equatorial  Undercurrent  do  not  show 
any  good  evidence  of  the  turbulent  subranges  required  for  the  Batchelor 
theory.  He  found  values  of  the  variance  of  temperature  fluctuations  apprec- 
iably lovv’er  than  observed  by  Williams  and  Gibson. 


Belyaev  ('t  al. (1975b)  have  reporteci  values  of  c for  llu-  Atlantic 
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Equatorial  Undercurrent.  They  observed  velocity  fluctuations  at  wavelengths 

loneer  than  2 cm  (kj~  3 cm  S and  estimated  values  of  e . The  values  are 

2 -3 

high,  ranging  from  about  0.  027  to  0.  34  cm  sec  for  horizontal  samples. 
There  is  no  significant  depth  variation  of  ^ over  the  depth  range  of  36  metres 
t(i  140  metres  which  they  sampled.  The  turbulence  is  very  intermittent,  and 
Ozniido\  (personal  communication)  has  indicated  that  the  values  of  t repor- 
ted represent  only  the  active  regions,  so  the  value  of  € averaged  over  all 
regions  may  be  a factor  of  ten  less. 

Measurements  of  temperature  and  microstructure  temperature  grad- 
ients by  Gregg  (1976)  at  155*W  on  the  equator  suggest  that  the  most  active 
turbulence  W'as  between  35  and  65  metres  depth,  which  is  above  the  velocity 
cor<’,  and  that  the  turbulence  intensity  at  the  core  is  much  less.  The  tur- 
bulence intensity  and  rate  of  dissipation  are  related;  hence  values  of  C 
where  Gregg's  data  were  collected  would  probably  be  relatively  high  above 
the  core  and  relatively  low  in  the  core.  The  sensors  employed  by  Gregg  were 
a thermistor  mounted  on  the  nose  of  a free-fall  instrument,  and  another  on 
an  outboard  wing  which  rotated  about  the  instrument.  Comparison  of  temp- 
erature spectra  from  the  two  thermistors  revealed  that  only  in  the  thermo- 
stad  is  there  a strong  suggestion  of  local  isotropy  at  microstructure  scales. 


I 


L 


33 


3 EXPERIMENTAL  APPARATUS  AND  PROCEDURE 

3.  1 Apparatus 

The  results  obtained  in  this  thesis  are  based  mainly  on  measure- 
ments by  a free-fall  microstructure  recording  instrument  during  1974. 

This  instrument  was  built  in  1971  and  1972  by  Dr.  T.R.  Osborn  and 
given  the  name  "Camel".  It  was  designed  to  measure  velocity  fluctua- 
tions in  a column  of  water  in  the  ocean,  as  well  as  temperature  and  con- 
ductivity variations.  Early  designs  have  been  described  by  Osborn  (1974) 
and  Osborn  and  Siddon  (1975). 

The  Camel  (shown  in  figure  6)  falls  through  the  water  with  sensors 
in  the  nose  to  sample  vertical  profiles.  The  fall  speed  is  40  to  50  cm/sec 
and  is  controlled  by  adjusting  the  trimming  weights  and  the  release  weights. 
The  brushes  at  the  upper  end  were  added  to  increase  form  drag  as  the  in- 
strument fell.  If  a device  such  as  a blunt  collar  were  installed  rather 
than  brushes,  the  large  eddies  shed  behind  it  could  cause  the  instrument  to 
wobble  as  it  fell,  and  contaminate  the  velocity  signal.  A fall  speed  of 
40  cm/sec  is  obtained  when  the  Camel  w'eighs  12  Newtons  in  seawater. 
Variation  of  density  over  the  range  encountered  did  not  cause  a significant 
change  in  terminal  fall  speed,  although  about  20  metres  are  required  to 
reach  this  speed. 

The  weight  release,  recovery  and  communications  systems  have 
been  described  by  Osborn  and  Siddon  (1975).  During  a free-fall,  signals 
are  converted  to  standard  IRIC  FM  frequencies,  and  transmitted  up  an  XW  I 
(an  expendable  wire  length)  to  a tape  recorder  on  board  the  work  ve.ssel. 

A Set  of  FM  discriminators  and  a chart  recorder  are  used  to  obtain  a 
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real  time  display  during  a free-fall.  In  1974  the  depth  range  of  the 
Camel  was  limited  to  450  metres  by  the  XWL,  but  a modification  of  the 
XWL  allowed  depths  of  800  m to  be  obtained  on  later  cruises. 

During  free-fall, eight  parameters  were  transmitted  to  the  ship: 
temperature  and  its  gradient,  conductivity  and  its  gradient,  gradients 
of  the  two  horizontal  velocity  components,  pressure  and  orientation  of 
the  Camel  relative  to  magnetic  north.  The  shear,  conductivity  and  tem- 
perature probes  were  located  as  shown  in  figure  6.  The  time  gradients 
were  obtained  electronically  within  the  Camel  before  transmission  to  the 
ship  to  improve  the  signal  to  noise  ratio  at  high  frequencies. 

The  method  of  obtaining  the  measurements  is  as  follows: 

V elocity 

The  main  probe  on  the  Camel  is  an  airfoil  velocity  probe,  used  to 
sense  horizontal  velocities.  Signals  from  the  probe  are  differentiated  to 
give  vertical  gradients  (assuming  Taylor's  hypothesis,  discussed  below) 
of  horizontal  velocities,  in  other  words,  the  shear.  Accordingly,  the 
probe  is  called  a shear  probe.  The  probe  was  developed  by  Siddon  (1965, 
1969)  for  studies  of  turbulence  in  wind  tunnels.  The  modification  for  use 
in  the  ocean  has  been  described  by  Osborn  and  Siddon  (1975)  and  the  use 
of  the  probe  for  studies  of  oceanic  turbulence  has  been  described  by 
Osborn  (1974).  The  calibration  and  testing  of  the  shear  probe  is  descri- 
bed briefly  in  the  next  few  pages,  and  more  detailed  information  is  given 


in  Appendix  A. 
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The  probe  is  of  an  axisymnoetric  shape,  as  shown  in  figure  7.  It 
is  mounted  at  the  front  of  the  Camel,  which  falls  nearly  vertically  through 
the  water.  The  expected  tilts  of  the  instrument  are  less  than  one  degree, 
occur  at  low  frequencies  (Osborn,  1974)  and  will  not  interfere  with  the 
measurements  of  fluctuating  velocities  which  are  found  in  the  dissipation 
range  at  higher  frequencies.  Water  flows  past  the  probe  at  a speed  ^ 
at  an  angle  of  attack  oc;  these  depend  on  the  fall  speed  of  the  Camel  V, 
and  the  horizontal  flow  u relative  to  the  Camel. 

The  sideways  force  on  the  probe  is  (Allen  and  Perkins,  1952): 

P = A Sin  2 OC  = smoc  cos  oc  3_  i 

= p AV  u 

where  ^ 

A is  the  effective  area  of  the  probe 

(5  is  the  density  of  the  fluid 

V is  the  fall  speed  of  the  Camel 

u is  the  horizontal  velocity  relative  to  the  Camel 
Al  is  (V^'  + u^)  (see  figure  7) 

The  derivation,  testing  and  limitations  of  this  formula  are  described  more 
fully  in  Appendix  A. 

The  probe  is  constructed  from  a soft  rubber  compound  and  contains 
tv\  o piezoceramic  bimorph  beams  mounted  perpendicular  to  each  other,  as 
illustrated  in  figure  7.  Much  of  the  transverse  force  on  the  rubber  mount- 
ing due  to  the  flow  is  opposed  by  these  beams  which  generate  x’oltages  pro- 


portional to  the  two  components  of  bending  moment. 
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Figure  7 


The  shear  probe. 
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To  calibrate  the  probe,  each  voltage  is  sent  to  a preamplifier, 
and  the  voltage  from  each  preamplifier  is 

Uci  - 2 Si^Vu.,  Wcz  ~ 2 Z 

where  Sj  is  the  sensitivity  of  the  probe  in  volts /(dyne/cm  ). 

For  the  measured  profiles  a differentiating  preamplifier  was  sub- 
stituted, a band  pass  amplifier  added,  and  the  output  voltage  became 

IT 


3.  2a 


,-l 


where  Kj  is  a circuit  constant  with  dimensions  of  (time)  as  explained 
in  Ajjpendix  C. 

The  probe  can  sense  two  horizontal,  perpendicular  components  of 
velocity  as  it  free-falls  through  the  ocean.  These  components  denoted  by 
U|  and  u^  , are  aligned  with  the  Camel  body  which  rotates  slowly  as  it 
free-falls.  A compass  was  installed  to  determine  the  rotation  rate.  For 
most  profiles  this  rate  was  below  1 rotation  per  minute.  The  signals  ha\-c 
not  been  transformed  back  to  north-south  and  east-west  components;  hence 
U|  and  u^  refer  to  the  output  of  the  first  and  second  channels  of  the  probe 
preamplifier,  and  have  an  arbitrary  orientation  to  the  large  scale  flow, 
hut  this  arbitrary  co-ordinate  system  does  not  affect  the  results  obtained. 

The  shear  can  be  obtained  with  the  use  of  Taylor's  hypothesis, 
\vhich  is  valid  in  this  case  (Osborn,  1974),  by  which  spatial  gradients  can 
be  derived  from  temporal  variations  sensed  by  the  probe  moving  relative 
to  the  fluid  at  a speed  V. 


^ _ -J_  I 
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The  equation  3.2  becomes,  for  the  velocity  components  uj  and 
— -L/o)  Ua  — - U~e»a  Kz 


3.  4 


where  Uoi,  Uoz  are  the  output  voltages  to  the  chart  recorders  on  board 
the  ship  (I  have  not  included  in  these  formulae  the  effects  of  the  filters, 
but  these  effects  are  noted  in  Appendix  C). 

A rubber  based  compound,  PRC-  1202-Q  medium,  was  used  in 
the  probes.  It  has  a hardness  on  the  Shore  A scale  of  about  40.  A typical 
\ alue  of  both  S|  and  was  5.0x10  ^ volts /(dynes /cm^ ).  Tests  showed 
this  sensitivity  changed  by  about  1%  per  degree  Celsius.  There  is  also 
a pyroelectric  effect  in  the  probes  by  which  a sudden  change  in  tempera- 
tvire  can  induce  a large  voltage  output  from  the  beams.  Because  the  sig- 
nal is  differentiated  before  transmission  to  the  ship,  it  is  the  rate  of 
change  of  voltage  which  is  important.  For  these  probes,  a step  change 
of  1“  C produces  a slow  voltage  change  of  100  millivolts /sec  which  corres- 
ponds to  a shear  of  0.  35  sec"^  at  a fall  speed  of  40  cm/sec.  However, 
the  output  signal  is  at  frequencies  low  compared  to  those  which  dominate 
contributions  to  the  velocity  derivatives  and  the  effect  was  reduced  by  a 
high  pass  filter  in  the  amplifier  (rated  -3db  at  \Hz)  (A  full  description  of 
the  electronics  for  the  shear  signals  is  given  in  Appendix  C).  Achlitional 
digital  filtering  was  performed  later  for  regions  of  low  turbulence-  inten- 
sity and  large  temperature  changes  to  further  reduce  noise  at  frequencies 
below  lllz,  and  so  prevent  leakage  into  the  frequency  batid  of  int(-rcst 


vs.  hen  the  Fourier  transforms  were  computed  (see  Apps-ndix  B). 


I. 


The  small  wavelength  response  is  determined  by  the  size  of  the 
probe.  -Siddon  (1969)  found  no  (i.  e.  ,< 0 . 5db)  attenuation  due  to  spatial 
averaging  of  an  earlier  model  of  the  probe  for  fluctuations  of  wavelengths 
longer  than  four  times  the  effective  length  of  the  probe.  Because  the 
pressure  forces  are  concentrated  near  the  tip,  the  diameter  is  chosen  as 
a measure  of  the  effective  length.  However,  probes  used  for  our  measu- 
renients  (diameter  4.  7 mm)  were  longer  than  the  probes  tested  by  Siddon 
relative  to  the  diameter,  and  the  effective  length  may  be  greater  than  the 
diameter.  A spectral  test  for  spatial  resolution  is  discussed  in  section 
4. 

To  determine  the  effect  of  high  pressure  upon  the  sensitivity  of 
tfie  probe,  the  Camel  was  launched  with  a bar  in  front  of  the  shear  probe 
to  generate  eddies.  The  Camel  descended  at  a uniform  rate,  so  the  tur- 
bulence behind  the  bar  would  have  been  uniform  with  depth  and  much 
larger  than  in  situ  turbulence  (the  gain  of  the  circuit  was  changed  to  avoid 
overloading).  No  change  in  output  signal  level  from  the  shear  probe  with 
depth  was  observed,  so  it  seems  that  the  sensitivity  is  independent  of 
de  pth . 

T empe  ratu  re 

Temperature  is  sensed  with  a Thermometries  microbead  thermis- 
tor which  is  nominally  0.013  cm  in  diameter,  coated  with  0.0016  cm  of 
F^aralene  C to  insulate  the  beads  anti  wire  leads  from  seawater.  Tests  in 
a flow'  of  12S  cm/sec  showed  the  response  was  reduced  to  half  power  at 


frequencies  of  15  Hz  (Lueck  et  al,  1976).  They  showed  the  response  was 
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mainly  determined  by  the  coating  of  Paralene  C which  is  a very  pood 
electrical  and  thermal  insulator.  Therefore,  at  slower  fall  speeds,  the 
response  time  should  not  be  much  different.  For  the  40  cm/sec  fall 
speed  of  the  Camel,  the  half  power  wavelength  is  about  3 cm. 

To  resolve  temperature  microstructure,  Gregg  (1973)  suggested 
that  wavelengths  of  1 cm  should  be  resolved.  This  resolution  could  be 
achieved  with  these  thermistors  if  the  response  function  were  known 
accurately,  since  it  is  a rather  simple  task  to  multiply  the  Fourier  trans- 
form of  the  temperature  gradient  signal  by  the  response  function  of  the 
niicrobead  to  obtain  the  spectral  form  of  the  actual  temperature  gradic'nt. 
This  procedure  has  been  used  by  various  researchers  to  resolve  the  t('m- 
perature  microstructure,  but  one  should  know  the  response  function  of 
each  individual  thermistor.  These  response  functions  were  not  availaljle 
for  our  measurements.  Also,  the  noise  level  of  the  temperature  gradient 
signal  (.002®C/cm)  may  interfere  with  signal  recovery. 

Another  factor,  often  neglected  with  these  thermistors,  is  fouling 
by  plankton.  When  using  the  response  function  of  a thermistor  to  recover 
the  temperature  gradient  between  3 and  1 cm,  a small  accumulation  of 
plankton  on  the  thermistor  can  change  the  response  function,  and  intro- 
fiuee  large  errors  in  the  calculated  gradient.  Fortimately,  it  has  no 
effect  on  the  average  temperature.  The  Equatorial  Atlantic  waters  are 
rc'latively  clear  of  plankton.  The  temperature  gradients  are  used  only 
in  a qualitative  way  in  this  thesis  and  plankton  fo\iling  prt'sents  no  proh- 
h'rn  for  this  purpose.  The  gradient  signals  havi'  not  been  corr('cted  for 
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the  rivgnal  attenuation  at  small  wavelengths. 

The  thermistors  were  calibrated  on  board  ship  prior  to  use.  Tem- 
peratures obtained  from  the  thermistors  used  on  the  Camel  matched  those 
of  the  CTD  within  0.  5°C  absolute,  and  their  precision  was  0.  2®C  over  a 
range  from  10"C  to  26®  C based  upon  comparisons  from  regions  of  small 
ti'mperature  change  with  depth. 

Conductivity 

A new  conductivity  probe  was  developed  during  1973,  based  upon 
till'  design  by  C.  S,  Cox  and  M.  C.  Cregg  at  Scripp's  Institute  of  Oceano- 
eraphy.  (Gregg  and  Cox,  1971).  It  provides  records  of  conducti\'ity 
gradient  in  the  microstructure  range,  but  the  conductivity  itself  showed 
a drift  which  was  depth  dependent.  At  the  equator,  conductivity  is  deter- 
miner! mainly  by  temperature,  and  the  two  signals  for  temperature  and 
< onductivity  grarlients  are  remarkably  similar.  Indeerl,  they  are  so  simi- 
lar that  only  the  temperature  gradient  output,  which  has  lower  noise 
le\els,  will  be  shown  in  the  discussion  of  the  results. 

The  ririft  of  output  voltage  of  the  conductivity  probe  with  fle]3th 
preventerl  an  accurate  determination  of  density  from  the  temperature' 
anrl  ( onrluct  ivity  signals.  Any  graphs  of  in  later  si'ctions  are  derix'i'd 
only  from  the  CTD  data  of  Dr.  Eli  Katz,. 

I 'ri's  su  re 

A vibrotron  installed  in  the  upper  end  cap  measure's  static  pressure' 
as  the  Camel  fre'e  falls.  The  frecpiency  of  a vibrating  wire  in  thi'  eihrei- 
t ron  \'.iri('s  about  a standee  rel  TRIG  fri'e|U('Mcy  as  the  amliient  ]iressure 
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alters  the  tension  of  the  wire.  A signal  of  this  frequency  is  transmitted 
to  the  ship,  giving  depth. 

The  fall  speed  is  computed  by  a computer  program  which  rejects 
spikes  in  the  digitized  record,  computes  a linear  fit  to  the  data  from 
successive  4 metre  sections  of  data,  and  calculates  the  fall  speed  ov'er 
the  section.  A detectable  fall  speed  change  (greater  than  3%),  likely  due 
to  a bubble  in  the  conductivity  bellows  was  observed  only  once  at  depths 
below  20  metres. 

The  vibrotrons  are  somewhat  temperature  sensitive,  with  the  cali- 
bration curves  shifting  for  changes  in  temperature.  Two  v'ibrotrons  were 
used  at  separate  times,  and  the  maximum  expected  error  in  average  fall 
speed  due  to  this  effect  is  1%. 

Direction 

A flux  gate  type  compass  installed  in  the  spring  of  1974  measured 
the  orientation  of  the  Camel  relative  to  magnetic  north,  and  hence  the 
orientation  of  the  individual  beams  in  the  shear  probe  relative  to  magnetic 
north.  The  compass  was  mounted  inside  the  Camel  pressure  case.  Out- 
put was  converted  to  an  FM  signal  and  transmitted  to  the  ship  with  tht' 
other  signals. 


1 


44 


.‘p*.*  '"L 


! 

i 


! 


I 


3.  2.  Procedure 

The  cruise  tracks  of  the  ATLANTIS  II  in  June  and  July  1974  are 
shown  in  figure  8.  The  peculiar  zig-zag  pattern  results  from  a combina- 
tion of  two  types  of  measurements.  The  short  paths  across  the  equator 
are  those  followed  when  an  instrumented  "fish"  was  towed  behind  the  .^hip. 
This  device  used  by  Dr.  E.  Katz  sampled  temperature,  conductivity,  and 
the  two  components  of  velocity  in  the  horizontal.  The  long  sections  at 
lO^W,  16“W,  22*V/,  28*W,  33“W  are  those  of  the  hydrographic  stations. 

The  Camel  was  used  on  these  sections,  along  with  a CTD,  Nansen  or 
Atlas  bottles,  and  current  meters.  The  data  to  be  shown  here  come  from 
the  sections  at  28®W  (profiles  20  and  22-25)  and  33°W  (profiles  28,  29)  as 
well  as  one  profile  at  Z4'’W  (profile  18).  For  most  Camel  launches,  a 
CTD  profile  and  current  meter  profile  were  taken  within  a few  hours  near 
the  same  spot  (within  a few  nautical  miles).  A current  meter  profile  was 
not  tak('n  at  24°W. 

The  Camel  samples  from  the  surface  to  between  300  and  400  metres. 
Reariings  near  the  surface  are  often  contaminated  by  oscillations  of  the 
Camel  caused  by  the  launching.  Low  noise  measurements  of  velocity  gra- 
rlient  iK'uin  between  8 to  20  metres  depth.  For  a Camel  launch,  tile  ship 
rlrifterl,  v^ith  starboarrl  sirle  up-wind.  The  Camel  was  hoisted  o\('r  the 
starljoarrl  side,  lowered  part  way  into  the  water,  and  released  froni  the 
launching  mechanism.  The  ship  drifted  during  the  free-  fall  of  the  Catiu'l 
to  acoid  cutting  the;  X\V  L on  the  ]jropellor.  A radio  on  the  Camel  signalled 
its  surfacing  but  the  direction  finder  did  not  work  wc'll  on  the-  steel -hu  1 led 
ship.  The  Caitiel  was  locatc-d  visually  by  looking  up-v\ind.  IL'iovery  of 
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the  Camel  was  achieved  by  snapping  a line  onto  the  recovery  hoops  and 
hoisting  the  Camel  on  board. 

There  were  twelve  free  falls  of  the  Camel  during  the  cruise. 

Those  with  reliable  shear  data  are  shown  in  section  4. 


I 


Figure  8 Cruise  83  of  ATLANTIS  II.  The  Camel  profiles  are  Indicated 

by  arrov's.  (from  Bruce  and  Katz,  1976) 
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Prof  i le 

Date 

Longitude 

Lat i tude 

18 

July  10 

24°0'W 

n°2'M 

20 

July  14 

28'’11'W 

2°o'N 

22 

July  15 

28°03'W 

n‘’37'N 

23 

July  16 

28°01'W 

0'18'S 

24 

July  16 

28‘’03'W 

n®17'S 

25 

July  17 

28‘’0’W 

1°21'S 

28 

July  23 

32‘’59'V 

0°2 

29 

July  23 

32°59'W 

P ° 2 ' N 

Table  I Camel  profiles  during  ATLANTIS  II  cruise 
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4 RESULTS 

4.  1 Qualitative  Aspects  of  the  Microstructure 

Chart  records  of  seven  of  the  observed  profiles  are  shown  in 
fiiiures  9 to  IS.  (A  current  meter  profile  was  not  taken  near  the  loca- 
tion of  profile  18,  and  this  profile  is  omitted,  although  dissipations  for 
profile  18  are  shown  later.  ) For  these  visual  displays  of  the  qualitative 
aspects,  the  traces  of  the  gradients  shown  have  not  been  corrected  for 
the  variations  of  the  sensitivities  of  the  probes  wdth  temperature,  and 
for  the  variations  in  fall  speed  of  the  Camel.  The  sensitivity  of  the  shear 
probe  changes  by  -7%over  the  range  of  temperature  found  through  the  pro- 
files, being  more  sensitive  in  warmer  waters  near  the  surface.  Thi'  ther- 
mistors are  more  sensitive  to  temperature  gradients  in  colder  waters  at 
rlepth,  changting  by  +4%.  Fall  speed  changes  below  ?0  metres  depth  are 
U'SB  than  3%,  except  for  profile  24  where  an  increase  of  7%  in  the  fall 
speed  was  found  between  the  depths  of  20  and  160  metres. 

The  shear  signals  have  been  band  passed  between  1 and  18  llz  (40 
anrl  2 cm)  to  reduce  both  low  and  high  frequency  noise  in  the  visual  dis- 
plays of  figures  9-17.  (These  filters  are  in  addition  to  the  filters  in  the 
elec  tronics,  described  in  Appendix  C.  ) The  temperature  gradients  are 
limited  by  the  response  of  the  thermistor  to  wavelengths  greater  than 
' about  1 cm,  and  have  no  low  frequency  (low  wavelength)  filtering.  The’ 

temperature  signals  have  been  low  passed  with  a 1 Hz  filter  to  reduce 
noise.  Moth  tempc’ ratvi  re  and  tempc’rature  graclient  signals  are  derivc'd 
frc)tii  thi'  same  probes.  Note  that  both  shear  and  tempe-rature  gradient 
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Profile  2.0  of  microstructurc  mcasuremenis  and  lemperatvir 
at  28“1I'W,  2®9'N  and  the  larpe  scale  current  profile  rnt-as- 
ured  by  J.  Rruce. 
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1 1 Profile  25  of  nue  rost  ruclure  mea  .su  rc'tnotit  s and  t empe  ralvi  r ( 

at  28®0'W,  1®21'S  and  the  large  scale  current  profih-  nieas- 
ured  by  .1.  Itruce. 
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sionals  include  most  but  not  all  of  their  respective  microstructure  ranges, 
with  the  small  scale  limitation  being  the  more  serious. 

All  the  profiles  except  number  20  were  taken  through  the  under- 
current. The  positions  of  the  Camel  profiles  are  marked  on  f igures  2 to 
4 to  give  an  indication  of  the  nature  of  the  undercurrent  at  each  profile. 

Also,  the  velocities  measured  by  Mr.  J.  Bruce  of  Woods  Hole  are  shown 
in  figures  9 to  15.  These  were  recorded  within  a few  hovirs  and  several 
kilometres  of  the  microstructure  profiles,  with  the  largest  difference  for 
profile  22  where  the  currents  were  measured  about  5 hours  before  and  5 
kilometres  away.  To  measure  these  currents,  a 6-10  minvite  average 
reading  of  a Braincon  model  252  meter  was  taken,  and  magnitudes  indica- 
tefl  are  relative  to  the  current  at  300  metres.  The  large  scale  velocities 
measured  by  the  current  nieter  were  not  sensed  by  the  slicar  probe,  which 
can  only  measure  velocity  gradients  at  mic rostruclu re  scales.  Gradients 
at  mic  rost  ructure  scales  are  usually  associated  with  turbulence-,  as  noted 
in  section  2.  Conversely,  the  microstructure  scale  shears  \sere  not 
sensed  by  the  current  meter. 

In  Table  II  the  difference  Richardson  numbers  (calculatcfl  using 
eciuation  2.  341  are  listed.  They  have  been  determini-d  from  the  current 
meter  readings  and  from  the  values  of  cr-j-  determined  by  Dr.  Katz  frorii 
the  CTl)  data.  These  values  should  be  treated  with  caution,  as  the  current 
nu-tc-r,  CTD  and  Camel  were-  ncjt  lowerecl  together,  but  were  (-ach  :-c'])arated 
in  time  and  space  from  the  other  by  sc-vc-ral  liours  ,i  nd  Kilometres,  Any 
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hil't  in  de[)th  c)f  the  n nch- r cu  r re  nt  betwi-eii  readings  of  current  and  density 
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coiilcl  causo  errors  in  the  \’alucs  of  Ri^  • 

riu'  profiles  Ihrouph  the  undercurrent  all  show  similar  fi-atures. 

I'he  most  intense  turbulence,  as  indicated  by  the  mean  square  microstrvu- 
"ir.-  \eU)city  erarlients,  is  usually  at  the  base  of  the  upper  mixed  layer, 
near  the  top  of  the  thermocline.  In  all  cases,  the  velocity  core  of  the 
ndercurri'iit  is  much  less  turbulent.  This  core  corrcsj>onds  rou}.dily  to 
ilie  r«.ee  ei  where  the  lari'c  scale  temperature  pradient  is  hif’hest.  ihdow 
ti'.e  ' elociiy  core  of  the  undercurrent,  the  turbulence  is  more  intcmse  than 
ai  t‘’e  core,  !)ut,  less  intense  than  the  turbulence  above  the;  core.  I'elow 
I'lis  tl-.ird  reeion,  the  profiles  are  nrarked  by  intermittent  turbulence  sepa- 
ralefl  liv  quieter  regions. 

Profile  dO  is  outside  the  undercurrent,  and  profile  d5  must  be  at  the 
extreme  southern  edpc',  as  the  peak  eastward  velocity  is  5 cm/scc. 

Profile  d5  shows  intense  turbulence  at  the  top  and  bottom  of  a very  sharp 
thermocline.  Profile  ?.0  has  its  most  intense  turbulence  in  the  the' rmocline. 

All  profiles  through  the  undercurrent  show  the  niost  active  turbulen- 
< e at  fh.'  to])  of  the  thermocline;  but  the  nature  of  the  turbulence  near  the 
(u  e,in  surface  varies  frotn  the  very  low  levels  found  in  profiles  23  and  24 
,il)o\i-  30  metres  to  the  quite  hiph  levels  found  in  profiles  28  and  20. 

Other  evidence'  showinp  the  most  active  lurbvih'nci'  to  h('  in  the 
upo’  r thertnocline  is  available.  For  example,  an  examination  of  the  tem- 
per,tiure  tide  rost  ructure  [X'ofih's  obtainerl  durinp  the  ATl.ANTkS  11  cruise 
b\  Mr.  I..  P.ilorh'au  usinij  sm.ill  free-falline  inst  r\iment  s , shows  that  all 
I !■  It  ])rofihs  lhro\u'h  the  u nde  r ( u r ri' nt  wi-ri'  .iitriilar  to  those  of  the  ('amel 
111  thi  upper  thermocline,  where  they  show  haree  I e m jX' ra  I n re  pr.idients  .it 
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niic  rost  rvicturc  scales.  Similar  features  were  found  in  the  Pacific  by 
('.reijc  (1976),  who  measured  temperature  in  the  Pacific  Equatorial  I'nder- 
I'urrent.  He  concluded  that  turbulence  was  most  int('nse  above  the  \-eloii- 
ty  (-ore  and  much  weaker  in  the  velocity  core.  It  appears  then,  that  the 
turbulence  levels  are  generally  highest  in  the  uppi-r  thermocline. 

The  large  scale  shears  have  values  as  large  as  .06  sc'c  ^ in  this 
rei'ion,  and  rms  turbulence  shears  are  up  to  . 3 sec  ^ averaged  over  10  m 
deiiths.  The  calculated  values  of  Ri^  in  the  upper  tliermocline,  shown  in 
Table  11  are  sometimes  below  0.25.  I’enci’  the  large  scale  shears  in  the 
upper  thermocline  of  the  undercurrent  may  be  sufficient  to  maintain  tur- 
bu  lence . 

In  all  profiles,  turbulence  shears  in  the  velocity  core  were  low. 
because  the  large  scale  density  gradients  are  high  here,  and  the  velocitv 
gradient  must  change  sign,  there  must  be  a maximum  in  the  dynamic 
stability.  The  magnitudes  of  the  turbulence  shears  increase  in  the  lower 
portion  of  the  thermocline,  below  the  velocity  core.  All  the  temperature 
microstructure  profiles  obtained  by  L.  Rilodeau  show  small  niicrostruc- 
ture  tempi'rature  gradients  in  the  middle  of  the  thermocline,  and  larger 
values  in  the  lovvt'r  thermocline.  The  values  of  Ri..  are  a maximum  in 
the  veloi  ity  core,  and  smaller  below  the  vudocity  ( Oi-e. 

I’lelovv  both  tlu’  Pai  ific  and  Atlantic  Ecpialorial  II  nrle  rcii  r rent  s 
rei-ion  of  nearly  uniform  temperatures  has  been  observaul,  with  1emi)era- 
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tureri  noar  IB^C.  This  feature,  referred  to  as  the  thermostad,  is  most 
evident  for  profiles  23  and  24  at  170  to  270  and  170  to  245  metres,  res- 
[)ec  ti\  ely.  The  value  of  Ri^  between  190  and  225  metres  is  0.3.  The 
accuracy  of  any  one  of  these  values  of  Ri^  is  not  sufficient  to  allow  one 
to  predict  whether  turbulence  may  be  present.  If  this  value  for  the 
Richardson  number  were  to  be  present  throughout  the  region,  then  shear 
instabilities  would  not  generate  turbulence.  However,  it  is  possible  that 
shears  are  concentrated  on  several  surfaces  within  the  thermostad,  and 
patches  of  turbulence  may  develop  from  shear  instabilities  as  noted  in 

section  2.6.  Below  the  thermostad  of  profiles  23  and  24,  intense  turbu-  j 

lence  is  found,  with  turbulence  shears  almost  as  large  as  those  in  the  ' 

upper  thcrmocline.  Other  profiles  do  not  show  thermostads  as  extensive 

as  observed  in  profiles  23  and  24,  nor  do  any  have  such  intense  turbulen-  j 

ce  Ijelow  200  n\ctres,  i 

Two  ])airs  of  profiles,  23,  24  and  28,  29  were  taken  consecutively.  i 

F’rofile  23  was  at  28®01'W,  0®18'S  at  1322  GMT  on  July  l6;  profile  24  was  ! 

at  2.8'’03'W,  0'’17'S  at  1435  GMT  as  determined  from  the  closest  satellite  i 

fixes.  The  turbulent  patch  size  and  intensity  of  the  two  profiles  are  sitni-  ; 

I 

j 

lar.  However,  the  depths  of  the  features  of  profile  24  are  1 0 to  30  metres  j 

i 

shallovvt'r  than  those  of  profile  23.  The  depths  of  features  of  profile  23  j 

ae  ree  to  within  a few  metres  with  those  recorded  by  the  CTD  meter  at 

’H^ni'W,  0®18'S  at  1219  GMT  to  1305  GMT.  There  may  be  some  distur- 
bance  or  wave  compressing  or  raising  the  featvires  of  the  water  column 
>,iii'pl<v!  by  these  proliles.  If  it  is  large  scale  and  long  jx'riod,  it  might 

I 
I 

L J 
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produce  quasi-steady  shears  which  could  maintain  a uniform  level  of 
turbulence  over  the  time  and  space  intervals  observed.  The  high  levels 
of  turbulence  observ€!d  in  the  thermostad  for  these  profiles  may  have  been 
enhanced  by  this  disturbance,  and  may  not  be  due  only  to  the  long  term 
• i\'erage  current. 

The  pair  of  profiles  ?.8  and  Z9  were  taken  near  32®59'W,  0“?,'X  at 
0917  CMT  and  1029  GMT  respectively  on  July  23.  Only  one  satellite'  fix 
was  available  for  this  pair,  but  1 believe  that  actual  positions  of  the  pro- 
files could  differ  by  several  kilometres.  For  these  profiles,  the  depths 
of  features  as  indicated  by  the  temp'-'rature  signals  are  similar.  Above 
the  c'orc,  the  turbulence  level  is  similar  for  these  profile  ..  In  the  lower 
thermocline,  the  region  o<'  turbulence  is  thicker  for  profile  29.  Below 
the  thermocline,  the  two  thin  patches  of  turbulence  at  280  metres  iind 
360  metres  observed  for  profile  28  have  decreased  in  intensity  for  pro- 
file 29.  The  noise  levels  of  the  gradient  outputs  near  18  1 z for  these  pro- 
files are  higher  than  found  on  other  profiles,  and  may  have  masked  some 
weak  turbulenct'. 

An  ('xpanded  c:hart  record  of  profile  23  is  shown  in  fieurcs  16a 
and  l6b.  These  expanded  gradient  profiles  show  the  relationship  betwu'en 
tlu?  velocity  and  temperature  mic  rosl  ructu  ri' , It  can  be  seen  that  \\her<’ 
temperature  microstructure  is  present,  and  both  ])ositi\e  and  negati\c 
microscale  temperature  g'radic’nts  are  obsi'rved,  there  is  always  velocity 
toic  rost  ructu  rc'  (and  henci'  turbulence).  I'urbulenie  requires  fhuluations 


of  all  three  (omponents  of  velocity  and  mixt's  a fluid.  Hence  a laree 
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Stall-  t I'lupe ratii rt-  ai-adient  will  be  broken  up  into  many  smaller  gradients, 
and  the  \t-rtical  turbulent  velocities  will  transport  colder  fluid  up,  warmier 
fluid  down  and  cause  the  positive  and  negative  temperature  gradients  obser- 


va-fl.  Where  temperature  microstructure  is  present,  but  only  positive  tern-  j 

perature  L’radientj  are  observed,  there  is  no  velocity  microstructure  ob- 
served. Examples  of  this  latter  type  of  micro-.!  ructure  are  found  at  69-70, 

'^  ! and  9b  metres. 

The  portion  of  the  record  near  69  metres  has  been  expanded  even 
l'■al■lher  in  figure  17.  I ere,  a large  scale  temperature  gradient  has  two 
ntaxima  which  extend  over  1.5  metres.  Many  sntall  positive  gradiento, 
s('pa  rated  by  about  4 cm,  are  found  within  this  large  scale  gradient.  The 
velocity  microstructure  here  is  below  the  noise  level.  In  the  patch  of 


iTtic  restructure  at  65  metres,  both  negative  gradients  and  some  turbulence 
ai-e  found.  This  region  from  65  to  80  metres  corresponds  to  the  core  of 
the  undercurrent,  where  one  expects  dynam.ic  -lability  to  be  the  greatest. 


I 

1 'the  interface  at  69  metres  also  marks  the  .'■alt  core  of  the  undercurrent. 


fh-low  this  point,  salinity  decreases  with  depth,  and  the  double  diff\ision 
regime  of  salt  fingering  is  possible. 

The  small  scale  positive  temperature  gradients  observed  at  69 
metres  may  be  the  remains  of  a mixing  event  for  which  the  turbulence  has 
decayed.  L.inden  and  Tvirner  (1975)  have  shown  that  internal  waves  can 
carry  away  frorn  a turbulent  region  a small  fraction  of  the  turbulent  ener- 
V if  there  is  a stably  stratified  region  above  or  below.  The  lurbuh-nce 
abov  e the  core  may  generate  such  waves  which  propagate  down  throiiuh 
the  core  anrl  '^etU'rate  instabilities  in  tne  core,  or  sharjicn  alreadv 
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rxistiriL'  temperature  gradients.  An  infinitely  sharp  temperature  gradient 
v.ill  spread  to  5 cm  after  about  an  hour  by  molecular  diffusion  (Osborn 
and  Cox,  1972).  This  length  corresponds  to  the  separation  of  these  small 

i 


gi-adient  spikes,  and  one  hour  may  be  taken  as  the  "age"  of  these  spikes, 
althoueh  it  may  be  an  ov'er-estimate  of  the  time  because  infinitely  slia  rp 
gradients  arc  not  expected  in  the  ocean. 

A second  possibility  is  that  these  small  gradients  may  be  genera- 
ted bv  an  interaction  between  turbulence  and  the  convection  caused  by 
salt  fingering.  I have  examined  the  expanded  microstructure  graphs  and 
have  found  that  temperature  gradients,  extending  over  more  than  50  cm, 
where  turbulence  is  below  the  noise  level,  have  these  small  positi\  e 
spikes  on  them.  Most  are  found  in  the  velocity  core  of  the  undercurrent, 
some  below  the  thermocline,  but  none  above  the  salinity  core.  Ik^nce 
they  arc  all  in  a region  where  salt  fingering  is  possible.  Unfortunately, 
little  is  known  about  the  interactions  between  turbulence  and  con\'ection 
set  up  by  salt  fingering. 

In  his  summary  of  measurements  in  the  Pacific  thermocline, 
Nasmyth  (1970)  noted  that  temperature  microstructure  was  often  found 
without  \elocity  microstructure  but  velocity  microstructure  (i.e.  , turbu- 
lence) was  always  accompanied  by  temperature  microstructure.  The 
flanu'l  profiles  show  regions  of  turbulence  where  the  large  scale  tempera- 
tiire  gradic-nt  is  small  and  temperature  microstructure  signals  lie  below 
tile  noisc'  level  (0.002  C^/cm)  of  the  circuit.  The’  ujspcr  mixed  layer  of 
profile  ’.5  abo\  i'  45  metres  shows  this  beha\  iour,  as  doc-s  a portion  of 
profile  ’,1  at  240  m,  in  the  thermostad.  In  almost  all  othi'r  regions  of 
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active  turbulence, temperature  microstr  icture  is  also  found. 

Because  turbulence  was  always  observed  in  the  thermocline,  above 
and  below  the  velocity  core,  it  may  be  that  the  seasonal  scale  velocity 
gradient  is  sufficiently  large  to  maintain  turbulence  in  this  region.  Even 
in  the  velocity  core,  the  "age"  of  the  temperature  gradients  is  less  lhan 
an  hour,  and  there  appears  to  be  some  process  present  to  stir  these 
waters  — either  shear  instabilities  or  double  diffusion  or  some  combina- 
tion of  these  two. 

4.  Z Spectra  of  Turbulence 

•Spectra  of  turbulence  were  computed  to  estimate  dissipation,  eva- 
luate the  degree  of  isotropy  and  provide  a useful  check  on  sensor  response. 
If  turbulence  is  locally  isotropic  at  wavenumbers  corresponding  to  the 
dissipation  range,  then  equation  2.  2Z  can  bo  used  to  calculate  viscous  dis- 
sipation from  a few  measurements.  If  turbulence  was  locally  isotropic 
over  a range  of  wavenumbers,  then  the  resolution  of  the  sensor  could  be 
determined  from  a comparison  of  the  known  spectral  shape  with  the  meas- 
ured. F'or  the  measurements  of  turbulence  using  the  Camel,  neither  the 
resolution  of  the  probe,  not  the  degree  of  isotropy  are  well  knov.n.  I ow- 
ever,  by  comparing  spectra  from  various  depths,  the  resolution  can  be 
determined  to  an  acceptable  degree  of  accuracy. 

There  arc  .several  profiles  of  microstructure  from  the  ATLANTIS 
II  cruise  which  show  regions  in  which  one  might  expect  sufficiently  locally 


isolrofiic  turbulence  to  lead  one  to  e.xpect  a universal  spectruni  (ner 


nicasurahlc  wavenumbers.  These  regions  are  the  upper  mixed  layers 
of  profiles  ?.5,  ?.8  and  29.  Here  the  stratification  is  low  and  the  mean 
sliear  is  high.  Difference  Richardson  numbers  approach  zero  for  sever- 
al 19  metre  thick  layers  in  these  profiles.  Of  the  two  scale  lengths 
'..I  'd  in  section  2.  S for  the  high  wavelength  limit  of  local  isotropy,  iliai 
of  fV.midoa  is  most  readily  calculated.  This  length  scale  is 

2.  32a 

For  thi.'  upper  layers,  this  scale  is  longer  than  50  cm  for  profile  25,  20 
cm  for  ]5r(tfile  28  and  15  cm  for  profile  29.  The  stratification  in  the 
uipier  layer  of  ]3rofile  25  is  so  small  that  measureable  temperature  and 
conductivity  gradients  at  microstructure  scales  have  disappeared  within 
the  upi>er  50  metres,  although  velocity  gradients  are  present  at  all  scales. 

The  power  spectra  of  velocity  components  for  those  regions  are 
shown  in  figures  18  and  19.  The  method  of  calculation  is  outlined  in 
Appendix  13.  The  spectra  are  from  three  sections  extending  through  about 
It)  metres,  and  one  through  abcjut  8 metres,  representing  8192  and  4096 
clata  jjoints  respectively,  and  are  band  averaged  as  explained  in  Appendix 
P).  They  have  been  fitted  to  the  universal  F2(k/k^)  curve  of  figure  5 
according  to  the  method  described  by  Stewart  and  Grant  (19621  which  is 
also  outlined  in  Appendix  B.  This  method  of  curve  fitting  yields  an  esti- 
mate of  th('  \alue  of  the  viscous  dissipation  £ , averaged  over  the  de])th 
CO  r fi- s [jond i ng  to  the  spectrum,  and  this  value  of  C is  denoted  as  Cu  to 
inrlicate  how  it  is  derived. 


Figure  18  shows  lh(>  fit  of  the  spectra  (normalized  using 


norr":al 'zed  sped  fa  vptsiis  log  (k/kg)  . 
n and  tiie  log  .iF  the  nifniilized  an 
nl  curve  as  described  in  Appendix  R. 


from  log  (k^/k  ) = 1.  3 to  -0.5,  which  encompasses  the  centre  of  the  dis- 
sipation spectrum  G2(k/k  ).  At  lower  wavenumbers  the  spectral  coeffi- 
cients fall  below  the  F2  curve,  due  to  a high  pass  filter  in  the  r hear  cir- 
cuit which  attenuates  the  signal  at  frequencies  less  than  the  half  power 

ko 

point  near  1 Hz.  At  values  of  log  (^)>-0.5  the  spectral  coefficients  fall 
below  the  universal  curve.  Beyond  the  last  point  shown,  noise  dominates 


the  spectra.  (The  noise  spectrum  will  be  shown  in  figure  ?.4.  ) The  out- 


put covers  a small  range  of  wavenumbers.  A broader  range  of  wavenum- 
bers would  be  desired  for  comparison  with  the  universal  curves. 

Figures  18  and  19  also  show  power  spectra  of  profiles  28  and  29 
where  more  active  turbulence  was  encountered.  Here  deviation  of  obser- 
ved  data  from  the  universal  curve  begins  at  values  of  log  ^ greater  than 

*^5 

-0.7,  just  beyond  the  peak  of  the  dissipation  spectrum  G^Ck/k^)  which  is 

k. 

nc.-a  r log  -0.8. 

In  the  power  spectra  of  figures  18  and  19  the  dcv'iatioii  of  tin.'  nu'as- 
ured  values  of  the  spectral  coefficients  from  the  universal  curve  is  about 
-1  db  {-201^)  at  11  Hz  (~4  cm),  and  -3  db  (-50"^,,)  at  15  Hz  (-3  cni). 

Ih-cause  this  behavior  of  the  power  spectra  coefficients  is  nearly  unitorm 
nvvr  a wirle  range  of  \ alucs  of  £ , I believe'  that  spatial  averaging  of  the 

probe  may  be  the  cause.  The  criterion  of  Siddon  (1969)  ir.  that  the  shear 
[ji-obe  can  resolve  to  wavelengths  as  small  as  four  tinu's  the  length  si  ale 

of  the  probe,  which  is  taken  to  be  the  diameter.  The  power  siK'(tr,»  of 
fii'iires  18  and  19  are  about  1 db  bi'low  the  uni\  ersal  cnrvi'  at  wavelengths 


of  eieht  times  the  protx'  dianieter.  Thi'  'hear  probes  used  in  this  stud\- 


were  proportionately  longer  than  those  of  oiddon,  and  it  may  be  that  these 
prohi'S  are  sensitive  to  eross-stream  flows  in  the  region  of  the  probe 
back  from  the  tip,  for  which  case  the  diameter  is  no  longer  the  appropri- 
ate length  'cale. 

Several  other  possibilities  have  been  examined  to  explain  the 

differences  between  the  observed  and  universal  spectra.  The  universal 

k 

curve  of  Nasmyth  was  chosen  for  computation  of  F2(T?s)  because  of  the 

lo\\  scatter  in  his  spectral  coefficients  over  almost  three  decades  of  uave- 

numbers,  the  extremely  large  Reynolds  number  of  the  turbulence  in  the 

tidal  channel  where  the  measurements  were  taken,  and  the  close  agree- 
k 

ment  of  his  F (Tc^  ) curve  with  those  obtained  by  others.  0\er  the  range  of 

non-dimensional  wavenumbers  used  for  comparison  in  figures  18  to  ?.  1 , 

k 

the  differences  in  shape  among  the  Nasmyth  curve  of  FCFg),  an  earlier 
cur\c  by  Crant  et  al.(19f>2),  and  later  results  by  Boston  (1970)  are  insign- 
ificant. Bioston  found  that  his  data  which  had  the  greatest  signal  to  noise 
ratio  agrc'ed  best  with  Nasmyth's  results. 

Tniversal  curves  might  depend  upon  the  degree  of  intermittence  of 
the  turbulence;  a dependence  most  noticeable  as  a change  in  curvature  of 
the  "knee"  of  the  spectra  (Pond,  I9f>5  and  Nasmyth,  1970  both  discuss 
this!.  I ha\a'  tested  my  spectra  with  a rough  test  used  by  Nasmyth  (1970) 
on  his  riata  anri  fo\ind  no  appreciable  difference  in  degree  of  intermittence 
betv.ct’n  the  two  sets  of  data.  If  indeecl  the  intermitlency  of  the  turbulence 
has  reduced  the  cuiwature  of  the  "knc'c"  of  the  cur\'e,  it  may  affect 


i 


N ismvih's  and  my  riata  in  a similar  way. 
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The  transfer  functions  of  all  the  filters  and  amplifiers  have  been 

checked  (see  Appendix  C)  and  their  high  frequency  roll-off  cannot  be  the 

cause  of  the  high  frequency  roll-off  of  my  spectra.  However,  a better 

agreement  can  be  found  for  these  spectra  if  the  actual  sensitivity  of  the 

shear  probes  is  greater  than  the  calibration  value,  A higher  sensitivity 

k 

would  give  closer  agreement  between  the  F2{l<g)  curve  and  calculated 

spectra  at  high  wavenumbers.  For  example,  the  separation  betw'een  the 

k 

calculated  spectra  and  the  curve  decreases  to  about  1 db  at  15  liz 

for  profile  29  if  the  sensitivity  is  increased  by  40%.  Such  an  increase  is 
greater  than  the  expected  error  in  the  calibration;  it  is  unlikely  that  this 
effect  can  explain  all  the  sudden  roll-off  of  the  spectra,  but  it  may  con- 
tribute some  apparent  roll-off. 

There  is  the  possibility  that  the  turbulence  does  not  follow  a uni- 
versal curve  at  the  wavelengths  sensed.  The  length  noted  at  the  beginn- 
inc  of  this  section  is  only  a guideline,  as  no  rigorous  tests  of  local  isotropy 
at  wavi'length  have  been  made.  There  arc  several  tests  for  local  iso- 
tropy wliich  can  be  applied  (equations  2.  18  and  2.  19),  which  test  for  the 
conflitions  upon  which  the  Kolmogoroff  hypotheses  are  based.  Of  these, 

I can  tt'St  only  the  relation 

For  the  spectra  :•  hov.n  in  figures  18  and  19  the  spectral  values  of  (p  ^ ^ 1 


and  <p22^^3^  appc'ar  to  b€‘  ne-arly  ec]ual.  Thi'  ratio 
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^ ■ ‘i-~rk:<pjk,)Tk, 

o 

determined  from  fitting  each  universal  curve  takes  values  of  1 . 1 , 1.1, 

1.0,  0.  7 for  figures  18  and  19.  These  values  are  as  close  to  1. 0 as  the 
errors  in  curve  fitting  will  allow.  This  test  gives  support  for  the  pre- 
sence of  local  isotropy,  but  does  not  prove  that  it  is  present.  The  spectra 

1 

do  follow  a universal  curve  at  wavelengths  greater  than  4 cm,  and  fall  i 

i 

below  the  universal  curve  at  shorter  wavelengths.  I cannot-  prove  that  the 
lalculated  spectra  are  of  locally  isotropic  turbulence  or  that  they  should 
follow  <i  universal  spectrum^ 

Such  spatial  averaging  of  the  probe  would  attenuate  fluctuations  at 

high  wavenumbers.  This  is  not  a problem  for  the  turbulence  found  in  the 

upper  layer  of  profile  25  (figure  18)  where  the  average  value  of  for  the 

-4  2 5 

(Pjj(k3)  spectra  is  2.7x10  cm  sec  . The  peak  of  the  dissipation 
spoctrutTi  is  at  k - . ^ which  has  a wavelength  of  10  cni. 

However,  the  larger  value  of  £u  determined  from  the  spectrum  of  profile 
29  from  35  to  51  m in  figure  19  is  7.  0x10" -’em  sec  , and  the  peak  of 
the  corresponding  dissipation  spectrum  is  at  4 cm.  This  peak  lies  close 
to  the  wavelengths  at  which  spatial  averaging  occurs,  and  velocity  fluc  tua- 
tions at  ::mall  wavelengths  (high  wavenumbers)  will  be  missed. 

In  figure  20,  a graph  of  C^(k/kg)  vs.  k/k^,  is  shown  for  the  three 

ri-gions  rc'presentc'cl  in  figures  18  and  19  which  have'  a large  range  in  Cu  . 

k 

[■'fniation  2.20  is  usi-d  to  derive  ('.^(17,)  from  the  measured  siic'ctra  of 

J 
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Figure  20  Normalized  dissipation  spectra  G^ik/k^)  vs.  f f r throe  of 

the  four  regions  of  figures  18  and  1°'. 
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(p  I 1(1x3)  and  Cp  22(^3)’  Cu  required  for  this  calculation 

is  determined  from  the  fit  of  the  <p^j(k3)  and  (p2Z^^3'>  spectra  to  the  uni- 
versal curve,  as  outlined  in  Appendix  B and  section  2.4.  The  value  of 
£u  is  not  influenced  by  the  spatial  averaging  of  the  probe,  and  viscosity 
chanues  are  small.  The  calculated  spectral  coefficients  of  G2(k/kg  )corr- 
e.spond  closely  to  the  values  determined  from  Nasmyth's  curve  (closed 
circles)  only  when  the  dissipation  rate  is  low',  as  found  in  profile  25.  At 
hicher  dissipations,  the  measured  spectral  values  beyond  the  peak  are 

smaller  than  expected  - -an  indication  of  the  degree  of  spatial  averaging 

- 3 2 - 3 

taking  place.  Note  that  even  at  a dissi.ation  rate  of  7.0x10  cm  sec 

the  i)robe  can  resolve  almost  to  the  peak  of  the  spectrum. 

In  figure  21,  I have  plotted  the  ratio  of  the  dissipation  obtained 

from  a fit  to  the  universal  curve  (1^  ) to  that  derived  from  the  measured 

\ elocity  gradients  (£  m ) , using  equation  B'14b.  The  absissa  is  the 

, 1 

Kolmoizoroff  microscale  = (£r»t^  )“*.  spectra  from  patches  of  turbu- 

lence below  the  thermocline,  as  well  as  the  upper  mixed  layer  were 
exaluated.  The  latter  are  found  much  closer  to  the  line  drawn,  which  is 
gi\fn  by 

Cu  - .fc58  V >2.8c-nn 

_ 4.2 

Cu  ^ 1.15  CtA  kfM  ''  2.8  cm'' 

I'he.se  form\ilae  wort;  used  to  correct  values  of  viscous  dissipation  Cn  . 
ihe  metliofl  of  calculating  £m  is  rliscussed  in  section  4.  3.  The  lower  limit 
of  ^ 1.15,  and  some  of  the  departure  of  data  points  from  the  fit  are 


lil  elv  due  to  thi'  method  of  calculating  . 


The  largi'st  departures  of  data 
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points  from  the  line  representing  the  correction  formula  are  50%  or  more, 
but  almost  all  values  of  ^ fall  within  ^30%  of  the  value  predicted  by  the 

‘ wH 

' formula. 

' Spectra  of  turbulence  in  regions  below  the  well  mixed  layer  are 

' shown  in  figures  22  and  23.  The  effect  of  stable  stratification  for  these 

profiles  may  be  considerable,  and  local  isotropy  may  not  be  found  even 

j 

at  the  w'avenumbers  of  the  dissipation  spectrum.  These  spectra  and 
' others  from  such  regions  show  more  variation  than  those  of  the  upper 

mixed  layer,  not  only  from  the  universal  curve,  but  also  from  each  other 
' and  this  variation  is  apparent  in  the  values  of  plotted  in  figure  23  par- 

ticularly  for  ksH’4  cm  Some  spectra  rise  up  at  low  wavenumbers, 
likely  due  to  the  pyroelectric  effect,  noted  in  section  3.  1.  The  greatest 
departure  of  these  spectra  from  the  universal  curve  is  at  low  wavenumbers, 
and  much  less  departure  is  found  for  log  "kg  > -0.  8 (i.  e.  , at  wavenumbers 
' beyond  the  dissipation  peak)  once  the  effect  of  spatial  averaging  is  account- 

' ed  for.  flinze  (1959,  p.  255)  shows  that  the  pressure  velocity  co-relation 

and  the  viscous  forces  are  more  efficient  at  higher  wavenumbers  in  equali- 
zing the  three  turbulent  velocity  components.  Hence,  the  tendency  toward 
isotropy  is  greater  at  high  wavenumbers  at  dissipation  scales. 


4.  3 Dissipation  Profiles 

For  the  digitized  data  of  profiles  during  the  cruise  of  the 
II  a fast-Fourier  transform  was  computed  for  every  consecutive 
points  (about  2 metres  depth).  The  input  signals  of  and  , 
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filtered  as  described  in  Appendix  C.  These  filters  consist  of  one  high 
pass  (-3db  at  1 Hz)  and  two  low  pass  filters  (each  is  -3db  at  about  60  Hz). 

In  addition  the  discriminators  are  -3db  near  80  Hz  and  60  Hz  for 
and  respectively,  and  roll  off  at  18db  per  octave  at  higher  frequencies. 

The  measured  dissipation  was  estimated  by  use  of  the  formula 


f - .ff—  z'c(j)c*(j)l 

^ ''ax3>'  Z ' ibOpS.V^/  L(/d2‘/)J  >j,  J B-1 


J=Jz 


4b 


-i'  2.  ' 2.60f>S,V^/  L(/d2‘/)J 

as  explained  in  Appendix  B.  The  values  of  and  are  the  low  and  high 

limits  of  the  transform  values  which  are  dominated  by  shear  signals  as 

explained  below.  A similar  expression  for  is  computed;  each  of 

■9x3 

these  estimates  of  dissipation  is  corrected  for  spatial  averaging  of  the 

probe  using  equation  4.  2,  and  the  average  of  these  two  corrected  values 

is  chosen  as  the  dissipation  of  the  two  meter  interval,  and  is  referred  to 

as  e . For  calculating  the  values  of  in  section  4.  2,  the  uncorrected 

Cm 

values  of  £„  are  used. 

The  values  C(j)  are  the  transform  values  returned  by  the  FFT. 

The  advantage  of  calculating  Gn  in  this  way  is  that  only  transform  values 

at  frequencies  where  signal  dominates  noise  are  included  in  the  summation 

(equation  B-14-b)  to  estimate  Cm  . To  determine  at  which  frequencies  the 

3 3 

signal  dominates,  the  spectra  of  log  (k^)  <p  and  log  (k3)  <p 

versus  log  f were  plotted  for  each  FFT  of  1024  data  points,  and  some  of 
these  spectra  are  shown  in  figure  24.  These  spectra  preserve  the  rela- 
tive' contributions  to  the  variance  of  and  (^*^^^3:3)  at  each  fre- 


cpiency,  and  were  employed  only  to  determine  those  relative  contributions; 


82 


hence  the  values  of  the  ordinate  axis  are  arbitrary. 

The  upper  spectrum  in  figure  24  is  typical  of  a region  of  very 

high  dissipation.  The  shear  signal  dominates  the  noise  at  all  frequencies 

up  to  about  30  Hz,  and  Fourier  transform  values  at  frequencies  up  to  30 

liz  were  summed,  although  values  at  frequencies  above  20  Hz  made  a very 

— 2 

small  contribution  to  the  total.  The  spectrum  below,  with  . 0021  cm 

sec  has  a peak  at  23  Hz  due  to  a vibration  of  the  probe,  and  coefficients 

— 2 

were  summed  only  up  to  20  Hz.  The  next  spectrum,  with  . 000025  cm 

. 3 

sec  , shows  a low  peak  in  the  shear  spectrum,  with  the  noise  dominating 
the  signal  at  frequencies  above  about  15  Hz,  and  coefficients  at  frequencies 
Ijeyond  15  Hz  were  not  included  in  the  summation.  Two  spectra  in  the 
lower  part  of  figure  24  show  regions  where  the  signal  was  weak  at  all 
frequencies.  At  some  depths,  where  the  temperature  change  was  large 
and  shear  signals  weak,  a high  pass  digital  filter  was  employed  (-3db  at 
about  1 Hz)  to  remove  low  frequency  energy  caused  by  the  pyroelectric 
effect,  which  may  have  "leaked”  to  higher  frequencies  in  the  FFT  compu- 
tation and  dominated  shear  signals.  The  dashed  line  in  the  spectrum  at 
the  bottom  of  figure  24  is  from  a record  w’hich  has  been  digitally  filtered. 
The  solid  line  is  a spectrum  of  the  record  before  digital  filtering.  There 
is  no  peak  in  either  spectrum  between  1 and  15  Hz,  and  no  dissipation 
coulfl  be  estimated.  It  was  fovind  that  for  spectra  where  no  peak  apj')i  • 
between  freciuencies  of  1 and  10  Hz  tlie  estimated  dissipation  has.  ■■ 
the  Fouri(’r  transform  values  between  1 and  10  Hz  (of  diei'.ilU  ■ 
data  if  necessary)  was  less  than  4x10  cm‘'sec  ' . ! ■ 
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Figure  24  Dissipation  spectra  of  various  Individual  2 metre  sections  of 
profile  23.  A digital  filter  removed  variance  at  low  frequen- 
cies in  the  dashed  line  spectrum.  The  units  on  the  vertical 
scale  are  arbitrary. 
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chosen  to  represent  the  background  noise  level. 

It  is  likely,  especially  for  regions  of  low  dissipation,  that  the  high 
pass  filters  removed  shear  signals  at  low  frequencies  because  spectra  in 
regions  of  low  dissipations  peak  at  lower  frequencies.  In  figure  21  it  can 
be  seen  that  the  values  of  Gj  exceed  by  about  15%  at  low  values  of  k,j^ 
and  hence  low  values  of  Gn  . The  estimated  value  Gu  is  not  disturbed  by 
filtering  at  frequencies  below  1 Hz,  whileGw  is,  so  the  difference  between 
these  two  may  be  due  in  part  to  the  filtering,  and  the  correction  formula 


(equation  4.  21  ) which  was  applied  to  each  individual  2 metre  estimate  of 
should  recover  some  of  the  lost  dissipation. 

The  high  frequency  limit  of  the  summation  of  the  C{j)C=''(j)  is 
generally  a function  of  the  rate  of  dissipation,  but  in  many  cases  it  was 
difficult  to  determine  where  the  shear  signal  ceased  to  dominate  spectra 
svich  as  the  middle  one  in  figure  24,  where  the  uncertainty  of  choosing  a 
high  frequency  limit  in  the  vicinity  of  10  Hz  results  in  an  uncertainty  of 
about  10%  in  and  adds  to  the  uncertainty  of  the  estimate  of  £ . How- 
ever, it  is  generally  a random  error  and  an  average  over  many  spectra 
may  reduce  the  uncertainty.  The  transition  between  signal  and  noise  in 
spectra  of  regions  of  higher  dissipation  is  sharper,  and  uncertainties  less. 

Portions  of  the  signals  of  arid  contain  spikes  which 

were  removed  from  the  digitized  record  before  the  FFT  computation  and 
replaced  with  zeros.  Often  these  spikes  were  very  narrow,  due  to  a brief 
dropout  in  the  signal  leading  into  the  discriminators.  Because  this  signal 
was  not  low  passed  filtered,  these  spikes  remained  narrow,  and  disturbed 
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little  of  the  surrounding  digitized  record.  Where  the  spikes  occupied  a 
significant  portion  of  the  record  for  either  of  the  or  signals, 

the  other  was  used  alone  to  estimate  "C-h  . If  both  signals  were  fouled 
extensively,  no  dissipation  was  estimated. 

The  plots  of  these  dissipations  are  shown  in  figures  25  to  30 
together  with  the  profiles  of  temperature,  salinity  and  sigma-T  provided 
by  Dr.  Katz  and  profiles  of  east-west  current  relative  to  the  current  at 
300  metres  provided  by  Mr.  Bruce.  The  scale  for  the  dissipation  is 
logarithmic.  The  magnitude  of  the  dissipation  in  the  2 metre  depth  inter- 
vals is  indicated  by  the  position  of  the  right  side  of  the  bar.  The  left  side 
62  3 

is  at  10  cm  sec  and  the  noise  level  of  the  system  corresponds  to  a 
dissipation  of  4x10  ^cm^sec  In  regions  where  no  bars  have  been 
drawn,  both  shear  signals  were  extensively  conta*'  .;ed  by  noise.  This 
usually  occurs  just  below  a region  of  large  temx  .re  change. 

The  estimates  of£,.,  rely  upon  the  sensitivity  of  the  shear  probes 
and  the  kinematic  viscosity  of  seawater,  both  of  which  change  with  temp- 
erature (the  pressure  effects  may  be  neglected).  To  account  for  these 
changes,  the  temperature  of  each  2 metre  spectrum,  determined  from  the 
thermistor  output,  was  used  to  correct  the  values  of  sensitivity  and  vis- 
cosity. 

The  viscosity  was  computed  from  the  formula  of  Miyake  and 
Koizumi  (1948).  Relative  changes  in  density  were  less  than  0.  5%  and  the 
average  density  of  all  profiles  was  used  for  the  estimates  of  Cm  . (The 
resulting  error  in  i,,  is  less  than  1%  and  has  been  neglected  . ) 
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(b)  Estimate  of  the  rate  of  viscous  dissipation  of  profile  18  at  24®0'W,  0*2'N, 
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Figure  29  (a)  Zonal  current  relative  to  300  metres,  salinity,  temperature  and  cr,  corresponding 

to  profile  25  (data  supplied  by  E.  Katz  and  J.  Bruce). 

(b)  Estimate  of  the  rate  of  viscous  dissipation  of  profile  25  at  28®0'\V,  1°21'S. 
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The  largest  factors  contributing  to  the  error  in  estimating  C are: 

(a)  The  errors  in  calculating  and  (^»  f derive  from  errors 

in  the  sensitivity  ( + 20%),  fall  speed  (+3%),  the  spectral  summation  (-10%) 

and  the  circuit  constants  K,  and  K-,  (^5%).  The  values  of  f and 

i ^ ^1X3'  3=C3'' 

vary  as  the  square  of  the  sensitivity,  the  fourth  power  of  the  fall  speed. 


the  first  power  of  the  spectral  summation,  and  the  square  of  the  circuit 

calibration,  and  assuming  these  errors  are  random  the  total  is  ‘ 

L(.h)^  +(.12/  - {.if  * CO"]  ""  = .95 

to  gi\'e  about  a 45%  error. 

(b)  The  errors  in  the  correction  formula  (equation  4.2)  based 
upon  the  scatter  of  the  points  appears  to  be  30%.  However,  the  value  of 
^ depends  linearly  upon  the  sensitivities  Sj  and  S2  and  the  circuit 
constants  and  K^.  A systematic  error  in  these  values  could  introduce 
additional  errors  into  the  correction  formula,  and  I have  allowed  a value 
of  50%)  for  the  total  errors.  The  correction  formula  boosts  the  value  of 

£ by  a factor  as  large  as  2.  5;  that  is,  it  accounts  for  as  much  as  60%  of 
the  final  estimate  of  the  dissipation,  and  a 50%o  error  in  the  formula  will 
give  a 30%)  error  in  the  dissipation. 

(c)  The  errors  dvie  to  the  assumption  of  isotropy  are  noted  in 
section  2.5,  and  may  lead  to  an  overestimate  of  almost  50%). 

A sum  of  the  errors  a+b+c  gives  a total  error  of  somewhat  more 

2 2 2 j 

than  a factor  of  2,  and  the  sum  (a^+b  +c  )^  is  somewhat  less  than  2 but 


it  is  unlikely  that  the  errors  would  combine  to  give  the  worst  case.  It 
is  also  unlikely  that  the  errors  are  random,  as  the  sensitivity  affects  the 
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errors  (a)  and  (b)  in  the  same  direction.  Hence  the  possible  error  in  the 
dissipations  plotted  in  figures  25  to  30  is  a factor  of  2. 

These  dissipation  profiles  put  in  quantitative  terms  the  profiles 
shown  in  fipures  9 to  15.  I have  discussed  the  upper  mixed  layer  and 
upper  thermocline  in  section  4.  1,  and  noted  that  it  is  always  turbulent 
above  the  undercurrent  core.  As  an  indication  of  the  relative  magnitudes 
of  the  dissipations,  1 have  calculated  average  dissipations  over  several 
regions,  and  compiled  them  in  Table  III.  Of  these  eight  profiles,  six 
are  at  or  reasonably  close  to  the  latitude  of  the  velocity  core.  Profile  20 
is  outside  the  undercurrent,  as  indicated  by  the  velocity  profile  in  figure 
26,  and  profile  25  was  located  in  a region  where  the  maximum  eastward 
flow  had  a magnitude  of  5 cm/sec.  Accordingly,  these  two  are  not  rep- 
resentative of  the  undercurrent,  but  they  are  of  interest. 

The  question  of  how  representative  these  six  profiles  are  arises. 
These  are  samples  through  the  undercurrent,  and  although  more  profiles 
would  be  desirable  the  variation  in  rates  of  dissipation  above  the  core 
from  one  profile  to  the  next  is  not  large  compared  to  the  depth  variation 
of  the  rate  of  dissipation  for  a given  profile.  For  example,  the  average 
rate  of  dissipation  in  the  region  of  intense  turbulence  above  the  core  found 
in  profile  18  is  3|  times  greater  than  found  in  the  same  region  of  profile 
23,  but  200  times  greater  than  typically  found  in  the  core.  I noted  in 
section  4.  1 that  this  pattern  of  turbulence  was  observed  above  and  in  the 
core  by  other  investigators  (Cregg,  1976  and  Bilodeau,  personal  comniu- 
nication).  This  turbulence  may  be  maintained  by  the  Reynolds  stress  vor’c- 
ing  on  the  local  mean  shear,  and  may  not  be  intermittently  generated. 
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Tf  one  measures  a short  term  current  profile  through  the  under- 
current, say  of  the  type  taken  by  Mr.  J.  Bruce  on  the  ATLANTIS  II,  it 
will  not  differ  appreciably  (i.  e.  , by  more  than  50%)  from  a current  ave- 
raged over  a month  or  so  at  the  same  spot.  Hence  the  large  scale  shear 
which  generates  the  turbulence  is  greatly  influenced  by  the  long  term 
average  current,  and  variations  in  the  intensity  of  turbulence  (and  also 
the  rate  of  dissipation)  may  not  be  large. 

Thus  these  rates  of  dissipation  in  Table  III  should  be  representative 
of  the  undercurrent.  This  situation  is  different  than  elsewhere  in  the 
oceans,  w'here  fluctuations  of  currents  are  usually  much  greater  than  the 
average. 

These  rates  of  dissipation  are  lower  than  reported  in  the  literature. 

2 ^ 

The  highest  average  rate  is  0.007  cm  sec”"’  for  profile  18,  and  is  low'er 
than  the  rate  of  0.08  cm  sec”  reported  by  Williams  and  Gibson  (1974)  in 
the  Pacific  Equatorial  Undercurrent.  Moreover,  these  investigators  report 
similar  dissipations  at  0®and  at  1®N,  whereas  the  values  in  Table  III  suggest 
lower  rates  of  dissipation  away  from  the  equator.  Values  of  £ reported  by 
Belyaev  et  al.  (1975b)  are  of  order  0.  1 cm  sec”  , and  allowing  that  their 
reported  values  for  the  Atlantic  Equatorial  Undercurrent  represent  regions 
of  intense  turbulence  which  may  be  a factor  of  ten  higher  in  dissipation  than 

2 - 3 

an  average  over  a large  region,  their  values  become  0.  01  cm  sec  , which 
arc  greater  by  a factor  of  3 than  the  average  of  my  estimates.  How'ever, 
their  values  show  no  significant  depth  variation  of  e betw'een  36  and  140 


i 

I 

j 


metres  and  the  results  noted  in  section  4.  1 and  in  Table  III  indicate  that 
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the  rates  of  dissipation  do  vary  by  several  factors  of  ten  over  these  depths. 

To  attempt  to  establish  the  expected  magnitude  of  t at  the  equator, 

I have  determined  the  estimated  values  of  rates  of  exchange  of  energy 
among  the  various  terms  in  the  energy  balance  for  the  average  and  fluctua- 
ting motions,  and  these  are  discussed  in  the  next  section. 

4.  4 Energy  Balance  in  the  Undercurrent 
(a)  Balance  for  the  Fluctuating  Motion 

Some  of  the  terms  in  the  balance  of  energy  for  the  average  flow 
must  be  determined  from  the  balance  for  the  turbulent  motion,  and  this 
latter  balance  will  be  discussed  now.  The  full  equation  is: 

^ _ - 1 ^ U-.  ^ g.q.^  - lTo,  ^U;  -c  - P S /s  2 11 

dt  <5Xi 

and  an  average  at  the  latitude  of  the  core  over  the  meandering  period  of 
16  days  is  desired  so  that  ^ aO  • Table  IV  shows  values  of  ^Di/  derived 

3X2 

from  the  current  meter  profiles  of  Bruce  and  Katz  (1976)  and  labelled  ave- 
rage '■hear,  and  the  values  of  £ from  the  Camel.  Neither  is  a 16  day 
average,  but  if  the  individual  measurements  for  the  four  profiles  are  ave- 
raged, they  may  represent  more  closely  a long  term  average  at  the  core. 
The  depths  of  features  of  profile  24  do  not  match  those  of  profile  23,  nor 
those  of  the  CTD,  and  so  profile  24  was  not  considered. 

For  each  depth  interval,  I have  assumed  that  £ , and  have 

neglected  the  contribution  of  u 2U  ^ 3U  ^ / 3 x ^ to  the  turbulent  energy  produc- 
tion. Ih'nce  -UjU^=  £/(SUj/5x^)  and  the  computed  value  is  liste 
Over  most  depths,  the  value  of  is  ni  vi  c h 


in  Table  IV. 
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Profile  22 


depth 

3Uj 

e 

-pu^U3 

Vv 

(m) 

(sec~l) 

^sec3^ 

(xl03) 

,cm2 

fc2) 

(xlO^) 

(£E^ 

sec 

20 

0 r\ 

-.022 

.16 

-7.3 

.3 

30 

40 

50 

57.5 
65 

72.5 
85 

112.5 

-.023 

4.2 

-180 

8 

-.017 

1.5 

-88 

5 

-.057 

.037 

-.065 

.01 

-.049 

,030 

-.061 

,01 

-.027 

.091 

-3.4 

.1 

-.018 

.020 

-1,1 

.06 

.026 

,24 

9,2 

.4 

Profile  23 


depth 

dUi 

e 

-PU3U3 

Vv 

(m) 

(sec"l) 

/cm2 

sec^ 

(xlO^) 

.cm2 

fe2) 

(xl03) 

sec 

20 

30 

40 

50 

60 

70 

77.5 

98 

148 

-.005 

.009 

-1.8 

.4 

-.002 

.17 

-85 

40 

-.012 

2,2 

-180 

15 

-.042 

3.8 

-90 

2 

-.031 

.20 

-6.5 

.2 

.023 

,015 

.65 

.03 

.00024 

.037 

154 

6.4x10 

.0088 

.093 

11 

1 

Profile  28 


depth 

e 

-PU3U3 

Vv 

8x3 

(m) 

(sec"’) 

.cm2 

sec^ 

.cm2 

sec2 

sec 

20 

30 

40 

50 

57.5 
65 

72.5 
90 

107.5 
125 
140 

(xlO^) 

(xl03) 

-.012 

1.9 

-160 

10 

-.008 

7.0 

-880 

100 

-.040 

8.3 

-210 

5 

-.020 

,40 

-20 

1 

-.010 

.47 

-47 

5 

-.027 

.059 

-2.2 

80 

-.0063 

.032 

-5.1 

1 

.031 

.012 

.39 

.01 

.020 

.73 

36 

2 

.013 

2.0 

150 

11 

Profile  29 

depth 

e 

-PU1U3 

V 

V 

(m) 

9x3 

(sec~^) 

2 

fe3) 

/Cm2 

sec2 

2 

(—  ) 
sec 

20 

30 

40 

50 

57.5 
65 

72.5 
90 

107.5 
125 
140 

(xl03) 

(xl03) 

-.012 

1.6 

-133 

10 

-.008 

3.4 

-425 

50 

-.040 

7.3 

-182 

5 

-.020 

2.6 

-130 

7 

-.010 

.30 

-30 

3 

-.027 

.027 

-1.0 

.001 

-.0063 

.67 

-106 

20 

.031 

.015 

.48 

.02 

.020 

.69 

35 

2 

.013 

1.4 

107 

8 

Table  IV 


Values  of  3Ui/3x3  from  data  of  Bruce  and  Katz(197fi)  and  c , 
and  derived  values  of  -U1U3  and  , the  vertical  eddy 
viscosity . 


smaller  than  <^'  , and  u^u  will  also  be  smaller  (assuming  the  stress 

^ 3 

varies  as  the  shear).  There  are  some  depths  where  is  comparable 

to  o)yi  , but  these  are  in  the  core  or  shallow  depths  away  from  the  region 
of  highest  dissipation,  and  will  not  greatly  influence  the  evaluation  of  terms 
in  equation  2.  11,  which  is  discussed  in  the  following  paragraphs.  The  ass- 
umption that  production  equals  dissipation  over  the  depth  interval  may  not 
be  strictly  true,  as  the  other  terms  may  be  significant.  However,  produc- 
tion and  dissipation  must  be  the  same  magnitude  and  use  of  this  approxima- 
tion allows  estimates  of  other  terms  which  allows  justification  of  the  assum- 
ption. 

. 2 

Above  the  core  a typical  derived  value  of  jujU^I  is  then  about  . 25cm 

'y 

licc^  in  the  region  of  highest  dissipation  and  this  is  the  square  of  the  fric- 
tion velocity  u^.  Monin  and  Yaglom  (1971,  p.  401)  estimate  for  a turbulent 

“2  ^ ^ 3 

shear  flow  in  a stratified  medium  that  q U3  and  p are  of  order  u ;;;  . 

? 2 

Smith  (1974)  has  found  u^q  . 05u,Uj  in  the  surface  layer  of  the  atmos- 

2 3 

phere  and  since  25u.,.  his  data  yield  n^q  s:  u .j^,  a result  consistent  with 

3 3 -3 

Monin  and  Yaglom's  statement.  I shall  take  2u,.,  - . 25  cm  sec  as  an  es- 

2 

timate  of  the  terms  u^q  and  p'u^/^  . 

Gradients  of  quantities  in  the  Xj  direction  will  be  much  smaller 
than  gradients  in  the  x direction,  and  may  be  neglected  (except  for  U^dq 

which  must  be  similar  to  x^).  For  a coefficient  of  correlation  of 

— 2 2 2 

0.  1 between  u ^ and  u^,  one  has  (uj  u^  )^  = 1 0 lujU^I  s?  2.  5 cm  / sec  and 

7 ■ ~ - 2 2 

(u  I )^si  (u^^")^  ^1,6  cm/sec,  (Actually  one  expects  (Uj  )>(n^  ),  but  for 
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A scale  height  may  be  taken  as  the  depth  from  a region  of  large 

€ to  one  where  £ is  low,  about  20  metres,  and  a mean  velocity  scale 

of  about  80  cm/sec  may  be  assumed;  and  choose  10*longitude  as  a zonal 

-3 

length  scale  and  1 day  as  a time  scale.  The  value  of  is  about  10  cm/ 
sec  and  will  be  discussed  in  the  next  section.  The  buoyancy  term  is 
assumed  about  0.  15  of  the  production  term  (Turner,  1973)  although  ex- 
perimental evidence  in  free  shear  flows  is  sparse  . (The  ratio  of  buoy- 
ancy to  production  terms  is  the  flux  Richardson  number,  discussed  in 
section  2.  5.  ) Equation  2.  11  is  now 


_ - 1 UjP'  . 

4 

IO‘'»  4 

80X4 

0.2‘b 

3x(0'^ 

'O* 

\o® 

2x  to^ 

Z *■  to"' 

I have 

chosen 

most  scales 

to  give 

an  upper  limit  to  the  terms 

compared  to  £ . For  example,  0.  1 is  probably  a low  correlation  coeffi- 
cient. In  the  atmospheric  surface  layer  a value  of  0.  3 is  typical.  The  ad- 
vective  terms  are  so  small  that  the  value  used  is  not  critical.  It  is  concei- 
vable that  the  buoyancy  term  could  be  up  to  . 2 of  the  production  term,  but 
this  value  is  not  significant  considering  the  uncertainty  of  the  values  e 
and  basic  balance  is  local  production  of  turbulence  equals 

local  dissipation. 

Below  the  core,  for  most  profiles,  £ is  smaller  by  a factor  of  3 
to  10,  but  other  terms  drop  correspondingly  and  the  balance  is  still  dissi- 
pation = production. 


(b)  Vertical  Eddy  Viscosity 


In  a turbulent  shear  flow,  the  Reynold's  stress  is  often  expressed 


(Hinze,  1959,  p.  466)  wherel^  is  the  vertical  eddy  viscosity.  Hence 

= _£ 4. 


I have  calculated  values  of  v>^  for  the  depths  between  the  current  meter 
readings  for  profiles  22,  23,  28  and  29  in  Table  IV.  The  value  of  »;■„  peaks 
at  between  8 and  100  in  the  high  dissipation  regions,  and  decreases  in 
magnitude  near  the  core  to  values  less  than  0.  1.  The  magnitude  increases 
somewhat  below  the  core  to  values  ranging  between  0.4  and  11  cm^sec"^. 

It  is  expected  that  the  value  of  vertical  eddy  viscosity  decreases  as  stra- 
tification increases  (Turner,  1973)  and  the  lower  values  of  in  the  core 
correspond  to  regions  of  largest  average  density  gradients.  Thus  it  does 
not  seem  possible  to  assign  a single  value  to  vJ,  for  the  undercurrent.  How- 
ever, in  the  region  of  large  dissipation  and  large  average  shear,  the  order 
of  magnitude  of  v),  is  10  cm'^sec  . Much  of  the  large  scatter  is  due  to 
the  different  sampling  times  and  durations  for  the  large  scale  currents 
and  the  microscale  shears. 

(c)  Balance  for  the  Average  Motion 

Again,  a time  average  over  the  meander  period  is  desired,  and 
meridional  terms  (i.  e.  , in  the  direction)  are  assumed  small  at  the 
equator.  In  this  case,  there  are  observations  of  several  components  of 


terms  in  the  energy  balance  averaged  over  the  GATE  Phases  of  about  20 
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days  each.  A compendium  of  measurements  of  anomalies  of  dynamic 

heights  along  the  equator  for  each  of  the  three  Phases  of  GATE  has  been 

prepared  by  Katz  et  al  (1976)  and  the  first  two  Phases  are  illustrated  in 

figure  31.  The  dynamic  heights  shown  correspond  to  zonal  pressure 

gradients±|f  , of  -3.2  and  -2.2x10  ^ dynes/gmfor  the  surface  and  50 

. 5 

decibar  depths  during  and  before  Phase  I,  and  -7.3  and  -5.3x10  dynes/ 
gm  during  Phase  II. 

The  change  in  slopes  between  Phases  I and  II  is  large,  and  is  main- 
ly based  upon  measurements  west  of  35’W  during  Phase  II,  not  available 
during  Phase  I.  It  can  be  seen  that  the  slope  of  equal  pressure  surfaces 
decreases  with  depth  through  the  undercurrent.  A diagram  of  the  change 
of  zonal  pressure  with  depth  has  been  prepared  by  Bubnov  et  al.(1975) 
and  is  shown  in  figure  32  along  with  velocity  profiles  from  moored  current 
metres . 

The  equation  for  the  balance  of  kinetic  energy  of  the  average  flow 
is 

^ P * L);  UiU^  -t-  uL'vUj 

<3t  P dX;  (5Xj  ^ 

We  can  try  the  scales  estimated  in  the  treatment  of  the  fluctuating 
flow  to  determine  the  magnitude  of  the  divergence  term  in  equation 

2.  7,  however  the  individual  components  in  this  term  are  also  those  in  the 
last  term  iTui  ^ , and  only  the  position  of  the  derivative  is  changed  for 
the  terms.  Scaling  will  not  show  which  term  is  larger. 

To  handle  the  term  equation  2.  7 may  be  integrated  over  a 


dx 


0 -20  -40  -60 


•40  -20  0 20  40  60  80 


Figure  32  Vertical  distribution  of  (a)  mean  east-west  pressure  gradient 
between  23°30'W  and  10°W  and  (b)  mean  profile  of  the  zonal 
current  at  23°30'W.  (Diling  1975,  figures  provided  by  V. Bub- 
nov) The  change  in  velocity  for  Phase  II  is  likely  due  to  a 
shift  of  the  undercurrent  to  the  south. 
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depth  range  in  the  direction.  An  integral  from  the  core  of  the  under- 
current at  = a (about  90  to  60  metres  depth)  to  the  depth  where  Uj  = 0 
at  Xj  = b (about  35  metres)  may  be  considered  as  one  region,  and  an 
integration  through  the  South  Equatorial  Current  from  = b to  the  surface 
at  Xj  = c may  be  taken  as  a second  region.  Because  there  is  a surface 
slope  the  value  of  c is  not  zero,  and  is  a function  of  x^.  Before  integrat- 
ing, some  comparisons  will  eliminate  small  terms. 

\^^  « L), 

« i^Di  ~ 

« <^lJiU>^  ~ (jLl?  <<  ^ Ui  U1U3 

5 3:, 

a}  ~ III  ~ ix^ 

The  eqviation  is  now 

f (jxs  - r - r ciXj  -f  lJ|  U.U3I  fu.Ui  <^^5 

f <5x,  J ^-x.,  J 5x3  ^ J 

' 4.fc 

The  magnitude  of  the  last  term  is  about  the  same  as  the  dissipation, 
as  shown  in  the  last  sub-section.  From  Table  III,  the  average  integrated 

3 

value  of  € from  20  metres  depth  to  the  core  is  somewhat  less  than  12  cm 
-3 

sec  and  only  profiles  28  and  29  have  large  dissipations  above  20  metres 
3 -3 

depth,  so  12  cm  sec  may  represent  a typical  value  from  the  surface  to 
core.  There  is  no  clear  tendency  for  this  dissipation  to  be  concentrated 
more  strongly  in  either  the  undercurrent  or  the  South  Equatorial  Current. 
To  evaluate  the  terms  in  equation  4.  6,  average  values  over  the 

region  of  integration  should  be  considered.  I can  assume  that  an  average 

2 7 2 7 2 2 

value  of  Q between  a and  b is  80‘'/4  cm  /sec  as  Q is  ^ at  the  core 

2 

(x^  = a)  and  0 at  the  level  of  no  zonal  motion  (x^  = b)  and  the  variation  in 
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Q between  = a and  = b is  reasonably  close  to  linear.  An  average 

value  for  the  thickness  of  the  upper  portion  of  the  undercurrent  is  40 

b ^ 

metres.  The  term  *^*3  is,  assuming  a change  of  flow  of  10%  per 

(90^  -€2*  j 3 + 3-3 

month  for  the  undercurrent  — tL fL x 4x10  ~ -1  cm  sec  . A 

•5  K io‘ 

significant  change  in  transport  of  the  undercurrent  has  not  been  reported 
for  the  three  Phases  of  GATE.  I have  assumed  that  a change  of  more  than 
10%  per  month  in  the  flow  of  the  undercurrent  would  have  been  observed. 

In  the  South  Equatorial  Current,  a larger  percentage  change  may  have 
taken  place,  but  the  velocities  were  somewhat  lower,  and  this  term  is  of 
similar  magnitude. 

The  term-fe^  dx_  may  be  determined  from  the  graphs  of  U,  and 
■'p  ax,  3 ' ^ 1 

^ in  figure  32,  or  from  the  anomalies  of  dynamic  height  in  figure  31 
plus  estimates  of  Uj  from  other  sources.  The  pressure  gradient  is  deri- 
ved from  the  anomaly  of  dynamic  height  AD{P)  by  assuming  the  reference 

level  is  indeed  on  a constant  geopotential  surface;  then 
. pcb^ 


-fU^P  dx5  = aAD(P)  dp 


4.  7 


Q 

where  A D(P)  is  in  dynamic  centimetres,  and  a similar  expression  is 
found  for  the  integral  from  x^  = b to  x^  = c.  From  figure  32,  it  can  be 
seen  that  the  integrand  changes  sign  at  the  level  of  no  motion.  The  esti- 
mated values  of  this  term  over  the  two  regions  are 

r'^Ll.  ^ ci-X,  ~ A ^ 

LlS  i cm  sec 


'b"?  3X 

.b 

! 

e 


^ ciX;  ~ +5 


3 -3 

cm  sec 


based  upon  a simple  integration  of  the  values  in  figure  32,  which  have  been 
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modified  in  the  following  way:  the  values  for  Uj  in  figure  32  are  derived 
from  averages  of  three  current  meters  at  32'S,  40'N  and  01 'N.  The  ave- 
rage current  at  the  equator  may  be  greater  than  the  indicated  value  of  70 
cm/sec  in  figure  32  by  about  10  cm/sec,  and  accordingly  a speed  of  80 
cm/sec  has  been  taken  for  scaling  the  core  velocity.  Also  the  speed  near 
the  surface  may  be  too  high.  The  current  meters  were  attached  to  surface 
floats,  which  may  have  increased  apparent  surface  speeds  by  wave  "pump- 
ins.’,".  A typical  value  for  surface  currents  measured  from  the  ATLANTIS 
II  is  about  -30  cm/sec  and  this  value  has  been  used. 

Note  that  the  slopes  of  the  0,  50  and  100  dbar  surfaces  in  figure 
51  indicated  by  those  ships  (KURCHATOV  and  PASSAT)  which  obtained  the 
data  of  figure  32  lie  between  the  slope  indicated  for  Phase  I and  that  of 
Phase  II,  and  thus  are  representative  of  the  average  pressure  gradients 
observed  by  all  ships  during  the  first  two  Phases  of  GATE.  The  observa- 
tions of  figure  32  have  been  used  for  analysis  because  they  show  the  varia- 
tion of  ^ , with  depth  in  greater  detail. 

The  two  divergence  terms  may  be  considered  together,  and  may 

aj-=a(x.) 


be  rewritten  as 

b 


f dU.aMas  = d r UiQ^dx^ 
J dx.  I 


tux.) 
b 


' au.) 


)a(Xi5  r UiGlM 

.JXi  L J 


( ’ dlXi  = rLl3Q^l 

> L -laixo 

aixil 
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I have  denoted  the  level  of  the  velocity  core  as  a(Xj),  because  the 
core  rises  toward  the  east  in  the  Atlantic.  Bruce  and  Katz  (1976)  observ- 
ed the  core  of  the  undercurrent  to  rise  by  abovit  30  metres  between  33*  W 
and  10*W,  and  this  figure  is  supported  by  other  data  (Duing,  1975).  There 


is  no  evidence  that  the  lev'el  of  no  zonal  velocity  (x^  = b)  is  a function  of 
X|.  There  is  also  evidence  that  the  vertical  velocity  a^  the  core  is  the 
vt'rtical  velocity  of  the  core;  in  other  words,  there  is  no  net  transport 
upward  through  the  core.  The  salinity  maximum  associated  with  the 
velocity  core  would  likely  decay  at  a more  rapid  rate  than  observed  if 
there  was  a significant  upwelling  velocity  through  the  core.  Also,  one 
would  expect  the  isotherms  to  ridge  or  trough  at  the  depth  of  the  core  if 
upwelling  or  downwelling  were  significant. 

In  this  case,  the  upward  slope  of  the  undercurrent  core  toward  the 
t>asl  is  the  ratio:  _ L)^  and  two  terms  in  equation  4.  8 cancel 

2 2 

The  value  of  U^Q  at  x^  = b is  zero  because  Q 2:  0 there.  Thus,  there 


remains  onlv  the  term 


'I  mtx.^ 


Two  effects  contribute  to  this  term:  because  the  core  of  the  under- 
current riics  and  the  level  of  no  motion  remains  reasonably  unifortn, 

there  will  be  a change  in  the  value  of  the  integral  along  the  Xj  direction; 

2 

secondly,  the  a\'erage  value  of  U^Q"^  may  change  in  the  Xj  direction.  To 

evahiate  the  first  effect,  assume  that  the  thickness  of  the  portion  of  tlie 

undercurrent  between  x^  = a and  b decreases  by  30  metres  from  33*\V  to 

> 2 

lO^W,  but  ihi'  average  energy  density  there  remains  the  same:  (80‘~)  cm",  . 

4 sec' 

The'  terni  becomes 


-jL  fU,Q^< 

djc,  J 


=z  - >0^  ) ^ 2 cm^: 


Zt  K to* 
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To  evaluate  the  effect  of  a change  in  flow  speed,  assume  that  Uj 
changes  by  10%  in  10°  of  latitude  between  3 3°W  and  lO'W.  Reported 
measurements  during  GATE  show  no  clear  trend  for  a change  in  U j along 
the  equator,  and  I am  assuming  that  a change  of  more  than  10%  per  10° 


longitude  would  have  been  sensed.  The  term  in  this  case  is 
b 


c/x 


ru.Q^jx,  * ± [f 

■ J 


+ 1 3 

= - 1 cm  o 


a(*i) 


lO 


sec 


where  the  average  values  of  Uj  over  the  limit  of  integration  are  taken  as 
88/v/2  and  80//?..  Of  course,  at  the  eastern  and  western  extremes  of  the 
undercurrent,  U|=0,  and  this  term  will  become  larger.  In  the  South 
Equatorial  Current  this  term  is  an  order  of  magnitude  smaller,  because 


the  value  of  Uj  is  much  less. 


The  term  -UjU^u^may  be  determined  from  the  average  velocity 


at  the  surface  (Uj~  -30cm/sec).  A yet  to  be  published  survey  by  Dr.  A. 
Bunker  of  Woods  Hole  of  atmospheric  winds  at  the  equator  gives  an  ave- 


rage value  of  the  zonal  wind  stress  - pu^u^^-O.S  dynes/cm  (Katz  et  al, 


-2 


1976)  and  so  -UjU^^-.  5 cm  sec  at  the  surface.  Hence  -U^u^u^^lb  cm' 


-3 


sec  at  the  surface.  Atthecore,  UjU^  = 0 should  hold  because  j xg 


0 and  changes  sign  through  the  core.  At  the  level  of  no  zonal 


velocity  UjU^u^^O,  because  U j = 0. 


Table  V shows  the  estimated  values  of  the  significant  terms  of 
t'ciuation  4.6.  These  values  may  be  considered  typical  between  10°W  and 
33'’W  for  the  currents  in  July  1974.  I have  portioned  tlie  terms  between 
the  .South  Equatorial  Current  and  the  Atlantic  Equatorial  Undercvirrent 
above  the  core.  The  pressure  gradient  varies  by  about  100%  from  Phase  1 


to  Phase  II  in  figure  31  and  estimates  of  €.  are  available  only  for  the 
western  portion  (24*W,  28*W  and  33“W),  and  may  be  in  error  by  a factor 
of  2 and  other  values  are  not  likely  to  be  more  accurate. 
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It  is  difficult  to  estimate  the  magnitude  of  the  meridional  terms. 

1 \ 

The  averaging  has  been  considered  at  the  equator,  where  the  terms  — ^ 
and  are  zero,  and  for  the  case  where  a meander  was  present,  but 

dXj 

not  growing  in  amplitude.  The  meridional  terms  3Ll>a.u^  and  are 

not  zero  at  the  equator,  but  could  be  small  compared  to  terms  such  as 
^ . Rather  than  try  to  evaluate  these  terms,  one  can  consider  the 

rate  of  energy  loss  from  the  average  flow  to  a meander  which  is  growing 
in  amplitude.  These  sixteen  day  meanders  were  the  most  energetic  of 
the  large  scale  fluctuations  observed  during  CATE,  and  the  rate  of  energy 
flow  to  a growing  meander  may  give  an  upper  limit  to  the  meridional  terms 
at  the  equator. 

As  the  meander  in  the  undercurrent  grows  and  sweeps  back  and 
forth  across  the  equator,  it  reduces  the  average  velocity  observed  at  the 
core  from  abovit  100  cm/sec  to  about  80  cm/sec.  K the  meander  period 

of  16  days  is  a time  scale  for  growth  of  the  meander,  then 

V ^ H<iO^  = -3  cm^/sec^ 

a J <)t  ' 3 

represents  the  rate  of  energy  loss  from  the  average  flow,  integrated  from 
th('  core  to  the  depth  of  no  motion.  The  magnitude  of  the  term  is  one 
half  of  the  integrated  zonal  pressure  gradient  term,  so  one  can  expect 
that  in  some  cases,  the  meridional  terms  may  not  be  negligible. 

A theoretical  study  of  these  meanders  by  Pliilander  (1976)  indicates 
that  the  Atlantic  Equatorial  Undercurrent  is  likely  stable  to  meanders,  hut 
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the  westward  flowing  South  Equatorial  Current  could  be  unstable,  and 
"perturbations  generated  by  the  instability  in  the  surface  layers  will 
cause  the  subsurface  undercurrent  to  meander".  The  energy  flow  for 
this  meander  is  from  the  average  flow  of  the  undercurrent  to  the  meander. 

Although  the  meridional  terms  may  be  important  in  the  energetics 
of  the  undercurrent,  the  region  I have  chosen  for  consideration  at  the  equ- 
ator between  the  core  and  the  surface  is  where  the  zonal  and  vertical  terms 
will  be  relatively  large,  and  the  meridional  terms  relatively  small,  and 
except  for  the  case  of  a growing  meander,  meridional  terms  should  be  un- 
important. Away  from  the  equator,  and  at  the  depth  of  the  core,  these 
meridional  terms  must  contribute  to  the  energy  balance,  and  I expect  that 
some  of  the  current  meters  moored  during  GATE  may  give  estimates  of 
these  terms. 

The  terms  shown  in  Table  V for  the  balance  of  energy  of  the  ave- 
rage motion  indicate  that  the  wind  stress  puts  energy  into  the  pressure 
field  found  in  the  .South  Equatorial  Current,  and  a portion  of  it  is  dissi- 
pated. In  the  Atlantic  Equatorial  Undercurrent  above  the  core,  the  basic 
balance  is  between  energy  removed  from  the  pressure  field,  adverted  by 
the  flow  and  lost  by  dissipation.  The  values  in  Table  V do  not  balance 
exactly,  but  do  agree  within  the  error  limits. 

For  the  portion  of  the  undercurrent  below  the  core,  the  rat(’  of 
dissipation  is  smaller  than  the  rate  found  above  the  core.  The  zonal 
I)ressure  gradient  is  also  smaller  liere,  and  puts  energy  into  both  thi‘ 


viscous  dissipation  and  the  adve-ctive  acceleration  tc'rms  (this  region 
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becomes  thicker  downstream  and  flows  faster  and  both  effects  take  energy 
from  the  pressure  gradient). 

The  general  agreement  of  the  terms  in  Table  V is  encouraging,  and 
increases  confidence  in  the  degree  to  which  the  eight  microstructure  pro- 
files are  representative  of  the  nature  of  the  undercurrent  at  the  equator. 

I believe  this  is  the  first  time  that  measurements  for  this  type  of 
energy  balance  have  been  taken,  and  the  undercurrents  are  the  only 
currents  for  which  a balance  as  simple  as  this  may  be  considered.  In 
other  current  systems,  the  balance  of  energy  is  not  as  simple,  due  to  the 
relatively  small  values  of  the  mean  velocities  compared  to  fluctuating 
velocities,  and  the  difficulty  in  measuring  the  relatively  small  cross-iso- 
bar flows.  However,  it  is  worthwhile  to  measure  the  rates  of  dissipation 
in  the  oceans,  and  for  special  regions  such  as  the  upper  mixed  layer  it 
may  be  that  terms  such  as  may  be  computed  from  the  estimates  of 

e.  , and  the  scales  at  which  mixing  takes  place  may  be  determined. 
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SUMMARY  AND  CONCLUSIONS 

During  the  GATE  project  in  1974,  undertaken  in  the  tropical 
Atlantic  Ocean,  a free-fall  instrument  was  employed  to  measure  vertical 
microstructure  gradients  of  temperature  and  horizontal  velocity.  These 
gradient  measurements  are  at  size  scales  where  the  turbulent  energy  is 
converted  to  heat  by  viscous  dissipation,  and  temperature  gradients  are 
smoothed  by  molecular  diffusion.  A portion  of  the  rate  of  viscous  dissi- 
pation can  be  measured  directly  by  a probe  on  the  instrument  and  the 
remaining  portion  is  estimated,  and  because  regions  where  the  rate  of 
viscous  dissipation  is  large  are  generally  regions  of  large  turbulent 
intensity,  one  can  compare  the  relative  turbulent  intensities  at  various 
depths . 

The  microstructure  measurements  at  the  equator  show  the  most 
intense  turbulence  is  at  the  boundary  between  the  upper  mixed  layer  and 
the  undercurrent,  where  the  average  shear  is  large.  Through  the 

core  the  turbulence  is  much  less  active  and  is  distributed  intermittently 
through  this  region.  Below  the  core  of  the  undercurrent,  at  the  base  of 
the  thermocline,  moderately  intense  turbulence  was  always  observed. 

Below  the  thermocline  the  turbulence  is  of  variable  intensity,  usually 
weak,  but  occasionally,  as  found  in  profiles  23  and  24,  some  turbulence  as 
active  as  that  found  above  the  undercurrent  was  observed. 

Spectra  of  turbulence  in  the  upper  mixed  layer  show  a consistent 
roll-off  from  the  expected  universal  Kolmogoroff  curve  at  frequencies 


above  10  Hz  (wavelengths  less  than  4 cm).  I have  interpreted  this  roll- 
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off  as  due  to  spatial  averaging  of  the  shear  probe,  and  have  compared 
rates  of  dissipation  €,j  estimated  from  the  fit  of  the  measured  spectra  to 
the  universal  curve  with  the  dis sipation  Cm  determined  directly  from  the 
measurements  of^',  and  , and  found  €u  was  consistently  greater  than 
. In  regions  of  very  intense  turbulence,  where  dissipations  are  high 
and  scales  of  turbulence  in  the  dissipation  range  extend  to  very  small 
wavelengths,  is  large,  up  to  2.  5 for  Kolmogoroff  wavenumbers  of 

10  cm  A correction  formula  has  been  determined  to  boost  the  estima- 
tes of  dissipation  to  account  for  this  spatial  averaging.  Thus,  the  assump- 
tion of  isotropy  is  invoked  to  predict  the  magnitude  of  the  unmeasured 
terms  in  the  expression  for  € , and  also  to  compensate  for  the  high  w'ave- 
number  portion  of  the  and  signals  which  was  missed  because  of  spa- 
tial  averaging.  The  expected  error  in  £ is  a factor  of  2. 

It  is  possible  to  estimate  terms  in  the  energy  balance  for  the 
fluctuating  motion  and  for  the  average  motion  from  these  microstructurc 
scale  measurements  and  measurements  of  large  scale  features  by  other 
investigators.  When  the  terms  in  the  energy  equation  for  the  fluctuating 
flow  are  scaled,  it  is  found  that  a balance  between  local  production  and 
local  dissipation  of  turbulent  energy  exists  where  turbulence  levels  are 
high.  Other  terms  in  the  equation  are  about  a factor  of  ten  smaller  in 
niaunitudc.  This  balance  allows  estimation  of  the  loss  to  turbulence  in 
the  average  energy  equation. 

Measurements  of  large  scale  features  during  GATE  (as  noted  by 
Doing  et  al,  1975)  revealed  a 16  day  -Z  day  meander  of  the  undercurrent. 


As  measurements  by  many  investigators  extend  over  an  entire  20  day  or 
so  Phase  of  GATE,  it  is  possible  to  average  some  readings  over  the 
meander  period.  Katz  et  al.(1976)  have  summarized  the  zonal  pressure 
gradient  data  for  each  of  the  three  Phases  of  GATE.  Bubnov  (Duing, 
1975)has  averaged  the  zonal  velocity  at  23.  5®W  and  the  pressure  gradient 
between  23.  5*and  10*W  over  each  of  the  three  Phases  of  GATE.  These 
measurements,  together  with  the  microstructure  measurements  allow 
estimates  of  the  zonal  and  vertical  terms  in  the  equation  for  the  kinetic 
energy  of  the  average  flow,  and  terms  are  evaluated  in  the  region  of  the 
equator  above  the  core  of  the  undercurrent  where  meridional  terms  are 
expected  to  be  relatively  small. 

When  the  equation  for  the  average  flow  is  integrated  between  the 
ocean  surface  and  the  depth  of  no  zonal  velocity,  it  is  found  that  the  rate 
of  kinetic  energy  gain  from  the  zonal  wind  stress  is  balanced  by  the  rate 
of  transfer  of  energy  into  the  pressure  field  and  the  rate  of  viscous  energy 
dissipation.  An  integration  from  the  core  to  the  depth  of  no  zonal  velocity 
shows  that  the  rate  of  increase  of  kinetic  energy  from  the  zonal  pressure 
gradient,  which  now  drives  the  current,  is  balanced  by  the  rate  of  energy 
loss  by  viscous  dissipation.  The  advection  of  kinetic  energy  by  the  flow 
is  probably  smaller  but  may  not  be  negligible. 
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Shear  Probe  Calibration 

For  a slender  body  of  revolution  in  an  inviscid  flow  of  speed  "Lf 
and  ansjle  of  attack  cc  , (assumed  small)  Allen  and  Perkins  (1952)  give 
the  cross  force  per  unit  length 

f - 5-2- 

f cross  force  per  unit  length 

f density  of  f 

rate  of  change  in  body  cross  sectional  area  with  longitudinal  dis- 
tance  along  the  body. 

The  authors  compare  this  theoretical  result  with  the  force  obser- 
vedona  model  of  a dirigible,  and  find  that  near  the  forward  tip  of  the  diri- 
gible, for  6*,  the  measurements  agree  well  with  theory,  and  viscous 
effects  are  small.  At  a = 1 2®  and  18®,  the  theory  underestimates  the  force 
slightly.  This  portion  of  the  dirigible  is  very  similar  to  the  tip  of  the 
shear  probe,  and  if  equation  A- 1 is  integrated  along  the  length  of  the  shear 
probe  from  the  tip  to  where  is  0.  then  the  total  cross  force  is 

r ~ sin  Zo<.  A-2 

(They  note  that  at  small  angles,  one  can  appro.ximate  sin  2a  ~ 2>x,  and  this 

approxitriation  has  been  made  by  Simpsoti  (l  972).  Tliere  is  a st-condary 

cross  force  which  varies  as  sin'^ot  and  is  felt  all  along  the  body,  not  just 

at  the  tip  where  0 and  this  force  has  been  neglected  for  small  \ ahies 

ax 


of  oc . 
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The  configuration  ofot.TJ  and  the  shear  probe  were  shown  in  fig- 
ure 7.  The  ability  of  the  shear  probe  to  measure  the  cross-stream  velo- 
city u is  evident  if  equation  A-2  is  rewritten  using  the  double  angle  rela- 
tion for  sin  2a. 

Fc  ~ C a A Z Sincx  cos  oc 

where  V and  u are  shown  in  figure  7.  It  can  be  seen  that  the  force  is 
linear  with  u. 

Another  approach,  described  by  Osborn  and  Siddon  (1975)  also 
yields  an  output  linear  with  u,  but  in  a slightly  different  way.  They  assume 

Tc  ~ rievMA(da\cx  a-4 

^doc  ! 

Where  is  the  change  in  lift  coefficient  with  angle  of  attack,  and  is 
doc 

assumed  constant  for  wavelengths  much  larger  than  the  effective  length  of 
the  sensor.  Equation  A-4  is  derived  from  an  extension  to  airfoil  theory  for 
a finite  width  airfoil.  At  small  angles,Ua  t\/u  and  this  expression  for 
approaches  the  form  of  equation  A-3. 

The  two  perpendicular  components  of  the  cross  stream  velocity 
are  sensed  by  orthogonal  biomorph  beams  within  the  probe  which  generate 
voltages  proportional  to  the  cross  force.  To  calibrate  the  probes,  these 
voltages  are  passed  to  identical  calibration  preamplifiers  whose  output 
voltages  are  recorded. 

Two  water  flow  devices  for  calibration  of  the  shear  probes  have 


bei-n  used.  One  operates  at  flow  speeds  of  100-150  cm/sec,  and  has  Ixu'n 
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described  by  Osborn  and  Siddon  (1975).  The  calibration  preamplifiers 
cannot  amplify  d.  c.  voltages,  so  a flow  is  directed  at  the  probe  by  a 
rotating  nozzle  inclined  at  5* to  the  axis  of  the  probe.  This  calibration  de- 
vice operates  poorly  at  flow  rates  below  100  cm/sec  and  rotation  rates 
less  than  16  Hz;  and  it  gives  no  check  of  linearity  of  response  with  the 
cross  stream  velocity. 

A second  calibrator,  shown  in  figure  A-1,  incorporates  major 
changes  to  allow  slower  flow  speeds  and  rotation  rates.  A jet  is  discharg- 
ed vertically  into  a tank  of  water,  with  the  tip  of  the  probe  mounted  jvist 
above  the  centre  of  the  outlet  at  an  angle  a.  Water  for  the  jet  passes 
through  a mesh  and  honeycomb  arrangement  to  reduce  the  scale  size  and 
intensity  of  turbulence.  The  honeycomb  also  reduces  large  scale  vortici- 
ty  which  would  produce  disastrous  effects  when  advected  through  the  red- 
ucer. A design  of  Smith  and  Wang  (1944)  was  used  for  the  reducer  to 
create  a uniform  flow  at  the  nozzle.  The  diameters  at  the  wide  and  nozzle 
ends  of  the  reducer  are  10  cm  and  2 cm  respectively,  to  give  an  area 
reduction  ratio  of  25. 

The  probe  is  tilted  successively  at  2|°  intervals  from  10®  left  to  10® 
right,  and  rotated  about  its  axis,  generating  sinusoidal  voltages  from  each 
bc'am  which  are  carried  to  the  preamplifiers  mounted  bcliind  the  probe. 

Slip  rings  bring  positive  and  negative  supply  voltages  from  the  preamp, 
atifl  tap  the  signals  from  the  two  channels  of  the  preamp.  A fifth  ring  is 
for  grounfl.  The  output  signals  are  sent  through  a low  pass  (-Sdb  a1  10  Ilz) 
unity  gain  filtcT  to  an  oscilloscope  to  check  for  any  bumps  or  spikes  in  the- 
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output  voltage.  The  filter  is  required  to  remove  high  frequency  noise 
generated  by  the  vibrations  of  the  slip  rings  and  of  the  electric  motor  used 
to  rotate  the  probe.  A rotation  rate  of  2.  5 Hz  was  used  for  calibrations. 

An  increased  rate  of  5 Kz  produced  no  change  of  output  voltage  magnitude. 

The  sinusoidal  output  voltages  from  the  filters  were  measured  with 
a true  rms  meter  (Disa  model  55  D35)  which  was  calibrated  with  a 2.  5 Hz 
sine  signal  whose  amplitude  was  determined  with  a well  calibrated  oscillo- 
scope. In  terms  of  the  rms  voltage,  the  peak  voltage  output  for  channel  1 
of  the  probe  and  calibration  preamplifier  is 
Ifc,  Lrr«>s  = 

= "s-pSiU^  SiT\«.cOS<X 

^ A-5 

= T ^syu 

where  Sj  is  the  sensitivity  of  channel  1 of  the  probe,  which  depends  not 
only  upon  the  nature  of  the  probe,  but  also  upon  the  gain  of  the  preampli- 
fier. A similar  expression  is  found  for  channel  2. 

On  board  ship,  a pump  was  used  to  drive  the  jet.  The  water  recir- 
culated through  the  system.  On  land,  an  overflowing  bucket  on  the  roof  of 
one  of  the  Oceanography  huts  at  UBC  provided  a uniform  pressure  head  to 
drive  (tie  jet.  The  use  of  this  calibrator  to  provide  accurate  values  of 
sensitivities  requires  the  consideration  of  several  factors. 

(a)  Gravity 

The  probe  tip  is  more  dense  than  water.  When  the  probe  is  rotated 


with  the  tank  full  of  water,  but  with  no  jet  flowing,  the  force  of  gravity 
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downward  on  the  probe  produces  a sinusoidal  voltage  at  the  output  ternain- 
als.  This  force  is  present  when  the  jet  flows  past  the  probe.  It  is  180® 
out  of  phase  with  the  force  of  the  jet.  The  ratio  of  the  voltages  from  the 
probe  caused  by  these  forces  at  40  cm/sec  flow  speed  in  the  calibrator  is 


Lj  is  the  lift  due  to  the  jet  (always  upward). 

Lp  is  the  gravitational  force  (always  downward). 

The  ratio  varies  with  V . This  large  effect  was  not  forseen  in  building 
the  calibrator.  With  the  jet  flowing,  the  net  lift  is  Lj  -Lg.  With  no  jet 
flowing,  but  the  probe  immersed  in  water,  the  net  lift  is  Lg.  The  corres- 
ponding output  rms  voltages  were  measured  at  various  angles  and  added 
together  to  produce  the  voltage  expected  from  Lj  alone.  This  procedure 
is  satisfactory  for  sinusoidal  signals,  but  any  noise  or  stray  signals  gene- 
rated would  be  added  doubly  by  this  method.  The  flow  speed  was  increased 
to  70-80  cm/sec,  where  |Lj/Lg|  s:  10  and  stray  noise  became  much  less 
troublesome.  Most  calibrations  were  at  these  flow  speeds. 

(b)  Boundary  layer  in  the  nozzle 

It  would  appear  at  first  that  the  flow  speed  through  the  nozzle  would 
be  the  volume  flow  rate  in  cm  /sec  divided  by  the  area  of  the  nozzle.  The 
flow  rate  can  be  measured  easily,  but  a boundary  layer  in  the  nozzle  may 
reduce  the  effective  cross-sectional  area  of  the  flow.  This  area  must  be 
known  to  caloilatc  the  flow  speed  U . The  sensitivity  varies  as  the  squari' 
of  this  speed,  and  the  measured  dissipation  varies  as  the  square  of 
the  sensitivity.  Thus  the  dissipation  depends  on  the  4th  power  of  effectivu' 
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cross-sectional  area  of  the  flow.  Small  boundary  layers  or  bubbles  on 
the  walls  of  the  nozzle  can  cause  large  errors  in  the  estimates  of  dissi- 
pation. 

Some  method  is  required  to  accurately  measure  the  speed.  The 

most  satisfactory  solution  found  is  to  measure  the  dynamic  pressure 
2 

in  the  jet  with  a manometer  and  pitot  tube.  This  system  has  the  ad- 
2 

vantage  of  measuring  V directly.  A large  bore  inclined  manometer  gave 
accurate  readings  of  the  height  of  the  water  column.  For  a pitot  tube,  a 
plastic  pipet  was  used,  with  part  of  the  tip  cut  away  to  leave  a 4 mm  end 
to  be  directed  into  the  flow.  The  dynamic  pressure  was  found  to  be  uni- 
form  across  the  core  of  the  jet.  Pressures  outside  the  jet  were  uniform 
also.  Abramovich  (1963)  states  that  observed  static  pressures  in  a sub- 
merged jet  such  as  observed  here  are  the  same  inside  and  outside  the  jet. 
Therefore,  the  difference  in  the  pressures  observed  is  the  dynamic  pres- 
sure of  the  flow  in  the  core  of  the  jet.  Figure  A-2  shows  the  flow  speed 

2 

for  various  volume  flow  rates.  The  value  of  U given  by  the  manometer 
should  be  accurate  to  better  than  1%  (Bradshaw,  1964).  The  repeatability 
of  results  was  within  2%  to  4%  and  should  be  the  error  in  measurements. 
There  is  no  reason  for  the  fit  to  be  linear,  but  the  line  drawn  does  fit  <^he 
points  well.  It  was  also  noted  that  the  flow  speed  increased  slowly  in 
time  d\iring  a steady  volume  flow  rate,  due  to  an  accumulation  of  bubbles 
I in  the  nozzle.  Flow  speeds  could  be  increased  by  5%  due  to  this  effect. 

During  the  ATLANTIS  II  cruise,  no  check  was  made  on  this  accumulation 
of  bubbles,  and  resulting  errors  in  sensitivity  can  be  up  to  10%.  For 
mf)re  recent  calibrations,  the  nozzle  has  been  brushed  clean  of  l>uhbles. 
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(c)  Turbulence 

The  large  scale  turbulence  was  reduced  by  using  packing  material 
in  the  bottom  of  the  calibrator,  and  aluminum  honeycomb  or  straws.  The 
honeycomb,  manufactured  by  Hexel  Corp.  was  3 mm  cell  diameter  and 
50  nim  deep.  It  was  mounted  9 cm  below  the  reducer.  The  Reynolds 
number  for  the  flow  through  each  cell  of  the  honeycomb  is 


Re  = = (o^ 


l_r  = 3 cm/s( 

d = 0.  3 cm 


U — 0.014  cm  / sec 

A-6 

I.aminar  flow  in  each  cell  can  be  expected  which  is  not  the  optimum  situa- 
tion for  the  generation  of  low  intensity  turbulence.  Turbulent  flow  in  each 
cell  would  be  preferred,  giving  lower  intensity  turbulence  beyond  the 
honeycomb.  For  laminar  cell  flow,  Lumley  and  McMahon  (1967)  predict, 
for  the  decay  of  turbulence  beyond  the  honeycomb 

^ ~ <e:£?d  a-7 

u'  turbulent  velocity 
U average  velocity 
d honeycomb  diameter 

X distance  downstream 

From  the  honeycomb  to  the  reducer  is  30  diameters  and  Ijp  should  be  0.001 
at  the  reducer.  Pankhurst  and  Holder  (195Z)  give  u'  ^ as  roughly  constant 
in  a reducer,  and  increases  by  (25)^,  pushing  turbulent  intensities  even 


lowt'  r . 


The  shear  proljc  itself  is  a good  detector  of  turbulence.  1 tric-d  a 


simple  test  to  determine  the  turbulence  level  in  the  flow  of  water.  With 
the  Jet  flowing  and  the  probe  inclined  at  10“  to  the  flow,  the  rms  voltage 
observed  with  the  probe  not  rotating  was  of  the  rms  voltage  observed 
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during  rotations. 

The  signal  from  the  rotating  probe  varies  as  U^'sin  PO**.  The  non- 

~2 

rotating  probe  will  not  detect  the  mean  cross  flow  because  signals  at  fre- 
quencies less  than  0.  1 Hz  are  filtered  by  the  preamplifier.  The  output 
signal  will  only  contain  the  higher  frequencies  of  cross  flows,  due  to  tvir- 

bulence,  which  vary  as  Vu'-T/U',  The  ratio  of  the  signals  is 
XJu' 

aj^sm20°  = 0-  0- 

from  which  /^/\j  ~0.  3%  which  indicates  low  levels  of  turbulence  at 
scales  and  frecjuencies  which  the  probe  and  filters  will  pass. 


(d)  Water  tunnel  blockage  effects 

In  the  ocean  the  probe  falls  through  the  ocean  at  speed  V . In  the 
calibration  tank  an  open  jet  of  w'ater  2 cm  across  is  directed  at  the  probe. 
The  deflection  of  streamlines  in  the  jet  is  not  the  same  as  in  the  ocean. 
The  jet  streamlines  spread  more  easily  around  the  probe,  reducing  the 
effective  dynamic  pressure.  If  the  probe  is  inclined  in  the  jet,  the  lift 
is  not  the  same  as  experienced  in  the  ocean.  Pankhurst  and  Holder  (195?) 


ha\e  described  the  following  corrections  for  these  effects  for  wind  tunnels. 

When  an  airfoil  is  introduced  into  a wind  tunnel,  the  solid  blockage 


at  the  mid-point  along  the  length  of  the  airfoil  is  conveniently  expressed 


Uf  - Ut  Cl 


A-H 
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U-p  tunnel  speed  far  from  the  model 

Up  tunnel  speed  at  the  midpoint  of  the  chord  of  the  model 

For  the  case  of  the  shear  probe,  one  wishes  to  find  the  blockage 
near  the  tip  of  the  probe,  where  ^ i-  0 (i.  e.  , where  the  probe  is  most 
sensitive  to  the  cross-flow).  One  possible  way  to  examine  the  blockage 
near  the  probe  tip  is  to  consider  only  the  rubber  portion  of  the  probe  as 
the  airfoil,  and  the  stainless  steel  tubing  behind  the  probe  as  a solid  wake. 
This  distinction  is  rather  arbitrary,  and  the  resulting  fineness  ratio  of  5 
is  higher  than  the  value  derived  if  only  the  sensitive  portion  of  the  probe 
is  considered  as  being  an  airfoil.  Pankhurst  and  Holder  have  cited 
results  of  Lock  and  Johansen  (1931  ) who  found  the  effect  of  wake  blockage 
small  when  compared  to  the  solid  blockage  of  an  airfoil.  In  the  absence 
of  better  information,  I assume  ..his  to  be  a valid  assumption  here,  and 
will  proceed  with  the  fineness  ratio  of  5. 

Pankhurst  and  Holder  (1952)  have  published  values  of  for  various 
flows.  For  an  open  axisymmetric  jet  blocked  by  an  axisymmetric  body. 


Es  is  given  by 


= ratio  of  cross-sectional  area  of  probe  to  cross-sectional 


area  of  nozzle. 


'J'  = constant  characteristic  of  nozzle  = -0.206  for  open  circular 


= c 


onstant  characteristic  of  airfoil.  This  is  chosen  as  5 which 


is  the  value  for  a Rankinc  Ovoid  of  fineness  ratio  5,  and 


4.0 
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Figure  A~3  Calibration  curve  of  the  shear  probe 


1 
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corresponds  to  the  assumed  fineness  ratio  of  the  probe  as 
noted  above. 

This  gives  fs  = C-o.20fc)(  5 )C  ^0.02. 

which  has  been  neglected.  A- 10 

The  second  tunnel  interference  effect  is  a change  in  lift.  Pankhurst 
and  Holder  (1952)  show  that  for  a finite  width  airfoil  in  an  open  jet,  the 
measured  lift  will  be  less  than  encountered  in  a flow  of  infinite  extent. 

The  change  in  lift  is  diminished  at  the  forward  and  trailing  ends  of  the 
airfoil,  so  for  the  case  of  the  airfoil  shear  probe  whose  lift  is  concentrated 
near  the  forward  tip,  this  effect  may  not  be  large.  The  corrections  given 
by  Pankhurst  and  Holder  apply  only  to  an  airfoil  whose  length  is  much 
smaller  than  the  jet  width,  and  do  not  apply  in  my  case  where  the  probe 
length  is  much  longer  than  the  jet  width. 

To  determine  the  magnitude  of  these  two  tunnel  interference  effects, 

I calibrated  a probe  using  two  different  nozzles,  of  diameter  1.92  and  3.00 
cm,  and  found,  for  similar  flow  speeds  that  outputs  with  the  larger  nozzle 
were  about  4%  to  5%  greater.  The  ratio  of  the  nozzle  areas  is  .41.  The 
blockage  effect  varies  inversely  as  the  three-halves  power  of  the  area  of 
the  nozzle  (equation  A-9)  and  the  lift  effect  varies  inversely  as  the  area 
(Pankhurst  and  Holder,  1952)  and  an  upper  limit  of  the  total  error  due  to 
these  two  effects  in  the  smaller  nozzle  is  found  by  assuming  the  sum  of 
these  two  errors  to  vary  inversely  as  the  nozzle  area.  In  such  a case, 
for  the  smaller  nozzle,  this  total  error  is  6%  to  8%  of  the  measured  lift. 

All  calibrations  were  done  with  the  smaller  nozzle,  and  no  correction  has 


been  made  for  this  error. 
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Probe  Sensitivities 

As  the  probe  rotates  in  the  jet,  the  lift  on  the  probe  remains  con- 
stant, but  the  lift  on  the  individual  beam  saries  sinusoidally.  The  peak 
output  soltage  is  for  channel  1: 

= 5,  {-p  Sin  2oc 

2 A-4 


In  practice,  the  rms  voltage  is  measured  and  the  peak  voltage 
derived.  To  find  S,  and  S,,  the  values  of  are  plotted  versus  sin  2ol 

as  in  figure  A-3.  The  flow  rate  is  measured  easily  by  trapping  the  over- 
flow from  the  tank.  The  speed  U"  is  determined  from  the  flow  rate  and 
the  graph  in  figure  A-Z.  The  output  voltage  is  linear  with  sin  Zoc  beyond 
oi-  Z.  5“  but  is  higher  at  small  angles.  This  deviation  is  due  to  non- 
perfect alignment  of  the  axis  of  the  probe  with  the  axis  of  the  jet,  such 
that  o<  = 0®  cannot  be  achieved.  The  relative  angle  of  the  jet  can  be  mea- 
sured accurately  (better  than  0.  Z°  ) but  the  absolute  angle,  as  shown  in 
figure  A-3  was  out  by  1®  . That  is,  the  reading  at  10®  to  the  right  was 
actually  at  1 1®  , and  the  reading  at  10°  to  the  left  was  actually  at  9°  , etc. 
One  would  e.xpect  that  because  sin  Zoc  is  not  a linear  function  ofocat  angles 
near  10°  that  an  error  would  result  from  mislabeling  the  angles  in  this  way; 
however,  the  error  is  small  ( < 1%). 

A larger  error  results  from  a tilt  of  the  entire  plane  of  inclination 
relative  to  the  jet  axis.  1 checked  this  angle  before  a calibration,  but 
could  not  be  sure  that  it  was  less  than  1®  . If  the  plane  of  inclination  is 
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figure  A-3 


Calibration  curve  of  the  shear  probe 


a nominal  angle  of  0“  . The  effect  of  this  tilt  is  diminished  at  larger  angles; 
the  true  angle  is  2%  larger  for  a nominal  angle  of  5“ , 1%  at  7j  and  0.  5%  at 
10“.  To  calculate  sensitivities,  most  weight  was  put  on  these  three  angles, 
and  the  expected  error  in  sensitivity  is  2%. 

The  largest  errors  result  from  uncertainties  in  the  flow  rate  in  ] 

the  calibrator  during  the  ATLANTIS  II  cruise,  and  the  tunnel  interference 
effects.  The  probes  were  recalibrated  upon  our  return  to  Vancouver,  but 
variations  of  about  10%  are  found  in  the  two  sets  of  calibrations.  It  may  be 
that  the  values  of  the  sensitivities  determined  during  the  ATLANTIS  II 
cruise  are  up  to  10%  too  high  because  of  bubbles  in  the  nozzle  during  the 
cruise,  and  6-8%  too  low  because  of  tunnel  interference  effects.  Other 
errors,  such  as  undertainties  in  the  flow  rate,  in  the  angle  of  tilt  and  in 
the  change  of  sensitivity  with  temperature  may  contribute  to  a lesser 
extent  to  the  overall  uncertainty  in  the  value  of  the  sensitivity.  I have 
allowed  that  the  overall  error  may  be  up  to  20%.  Errors  in  the  probe  sen- 
sitivities are  smaller  for  recent  calibrations. 


APPENDIX  B 


Data  Analysis 

The  analogue  signals  were  transmitted  from  the  Camel  up  a XBT 
wire  as  FM  multiplexed  signals,  and  recorded  on  a Hewlett  Packard  3960 
Instrumentation  Recorder.  A 14.  5 kHz  signal  was  also  recorded.  The 
signals  were  demodulated  and  demultiplexed  by  a set  of  Sonex  discrimina- 
tors, the  14.  5 kHz  signal  being  used  for  tape  speed  compensation.  The 
shear  signals  were  carried  on  the  2.3  and  3.0  kHz  channels.  These  ana- 
logue signals  were  digitized  at  200  times  per  second  onto  9-track  magne- 
tic tape  by  a 1 0 bit  converter  on  a PDP-12  computer.  Analysis  was  done 
on  an  IBM  370  computer  at  the  University  of  British  Columbia.  All  eight 
signals  from  the  Camel  were  digitized  in  this  manner.  The  horizontal 

velocities  u^  and  u^  were  sensed  by  the  probe  sampling  in  the  x^  direction. 

2 

The  corresponding  spectra  measured  are  ^jj(k^),  ^jj(k^) 

a n d ( k (p  2 2 ^ ^ 3 ^ * 


Dissipation  Calculations 

The  two  metre  averages  of  the  viscous  dissipation  were  derived 
from  fast  Fourier  transforms  (FFT)  of  1024  digitized  points.  Each  FFT 
gave  513  complex  transform  values,  and  the  square  of  the  moduli  of  the 
last  512  were  averaged  in  blocks  of  4 and  the  cumulative  svims  of  the 
first  45  blocks  (for  frequencies  up  to  35  Hz)  were  printed  out.  The  first 
transform  value  is  the  mean,  and  was  always  zero;  (the  average  was  com- 
puted and  subtracted  from  each  digitized  record  of  1024  points  before  the 
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I 

k 


FFT  computation). 

The  measured  signals,  transmitted  to  the  ship  and  recorded  are 


of  , <3uj  . The  values  are  related  to  the  output  voltages  by 
5-t 


_ 1,60  lToi 
cDt  fS.v 


B-1 


= 1.75  u'oz 


B-2 


as  shown  in  Appendix  C and,  by  use  of  Taylor's  hypothesis,  valid  in  this 
case,  one  has 


_ -1.80 


LToi 


(oS.yi 


B-3 


olUz  ==  -1.75  cToa 


B-4 


The  effects  of  the  high  and  low  pass  filters  have  been  neglected 


here,  but  they  are  discussed  in  Appendix  C.  The  constants  1.80  and  1.75 

-1 


in  equations  B-1  to  B-4  have  dimensions  of  seconds  , and  are  referred 
to  in  section  3.  1 as  2Kj  and  ^^2,' 

When  digitized,  one  volt  of  input  is  given  a value  of  250  bits;  thus, 
when  converting  the  digitized  output  to  volts,  one  must  multiply  by  (1  /250). 
Hence 


( = - I BP  An 

^ I S 250^5, 


50  (>S,' 


B-5 


where  d|  is  the  i -th  digital  value,  and  I ranges  from  0 to  1 023.  The  FFT 
floes  not  give  the  Fourier  sine  and  cosine  coefficients  directly,  but  the 


following  discussion,  from  Lee  (1974)  shows  how  the  discrete  Fourier 


transform  values  (i.  e.  , the  values  given  by  the  FFT)  are  related  to  the 
Fourier  series  sine  and  cosine  coefficients. 

Given  N samples  of  real  data  d (where  =0,  1 , . . , . N-  1 ) taken  at 
equally  spaced  intervals  At=T/N  (where  T is  the  period),  a corresponding 
Fourier  series  is 


■ ! C 

cl(-t)  = a.  +^(ajCos2i|H  bjsm 


If  t = 0 at  do,  then  for  each  d ^ we  have  t_=_^and 

T N 

din  ^ cto  + E (a, + b.s>h 

2 j=l  -v  N 


th°  FFT  of  the  d^  if  C(j),  and 


N ,-o 


It  is  possible  to  show  that  when  the  d(^)  are  real 

Md(J?)  = Co  + 2.|l|>e(c(i))c‘>s  - Im(CO^)s'ri^j-2  + (-0^C(%') 

j='  ^ R-9 

where  Re  is  the  real  part  and  Im  is  the  imaginary  part.  This  immediately 
gives  the  following  results; 

Co  = Qo  =0  in  our  case  because  the  average  is  zero 
Re(CC))  = ^Qj  Im(CO))=-^bj 

C(^"a)  = n-i 


n-io 


It  can  be  shown  that  the  variance  for  d is 
V-/  i-i 


M x-o  c-  j.i  -> 

- i[(j)'tc/;  .c4i 

That  is.  the  sums  of  the  squares  of  the  Fourier  series  coefficients 
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divided  by  2 equals  the  mean  square  value  of  d.  For  example,  if  the  time 
series  d(^)  is 


Ail)  = D cos 


2rr^ 


A/ 


B-  1 2 


^ z 

The  average  value  d will  be  D /2.  The  FFT  of  the  series  in  equation 

B-11  will  return  only  one  non-zero  coefficient:  C(l)  = ND,  and  the  aver- 

2 

age  of  the  power  series  will  be 


B-13 


because  C(l)  in  this  case  is  real. 
Accordingly 


= ' 5 v( 

^ 1.80  \^r 

a ' 

.250f.S.vV  L 

z 


p -I  B-14a 

and  a similar  expression  holds  for  In  practice  there  is  noise, 

particularly  at  high  frequencies  and  one  cannot  use  the  full  sum.  One  can 
however  examine  the  spectrum  to  determine  the  upper  and  lower  limits  of 
the  sum.  The  last  term  in  the  summation  of  B-14a  was  always  dominated 
l)y  noise,  and  if  the  upper  and  lower  limits  of  the  summation  for  the  shear 
spc'ctrum  where  noise  is  relatively  small  are  and  then  equation  B-14a 
b(.' comes 


^ 'asops.vv  L'"  J=J,  j B-14b 

It  is  the  separation  of  noise  and  signal  in  frequency  space  which  is  the 
advantage  of  calculating  i.  in  this  way. 


141 


Spectra  of  and 

The  formulae  for  are  given  below  and  the  formulae  for  <^tare 
the  same  except  that  the  constant  1.80  is  replaced  by  1.75.  To  determine 
the  relative  contributions  of  signal  and  noise  to  the  values  of  C(j)C*(i)  the 
spectral  density  functions  are  plotted  for  each  FFT  of  1 0Z4  data  points. 
The  spectral  density  function  of  ( ) is  (k^)  (p  ^.s  given  in  section 

2.  4. 


£ = 


B-  1 5 


If  C(j)  is  the  jth  transform  value  of  the  signal,  then  the  spectral 

(93t3 


density  function  is 

(k,J<p„j(h)  j=  >.2.  f . 

-(LS2. f r->-  cti>c’(j))  j*" 

^2/v2  ate,  / 2 


B-16 


The  is  in  the  denominator  because  each  Fourier  series  coefficient 
must  be  divided  by  the  portion  of  wavenumber  space  which  it  represents, 
= ^TT^f  _ 


Ii3i  = ^ 


VMAt 


B-r 


B-18 


The  spectra  of  figure  24  are  log  (k  ) (p  (k  ) versus  log  k,,  and  sub.sti- 

3j  ^ 1 Ij  3 ^ 

tuting  equations  B-17  and  B-18  into  B-l6. 

( (i,)  - 5-' 


I J C(s)c*(j))  , = ^ 

'ZSO^S.V'V  V2v^  / ^ 


B-  Ih 
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A plot  of  (k^.)^<p  versus  log  may  be  integrated  to  give  the  vari- 


ance of  the  signal: 


2 3 J (P>i  (ki)  ci  ^ ki  = /Yk) 

o 

- / (kf^Jkl)cJk, 


B-20 


3 

Such  a plot  is  variance  preserving.  In  my  case,  a plot  log  (k3j)  <pjj.(k^) 
is  used,  which  is  not  variance  preserving,  but  one  can  see  the  relative 
contributions  to  the  variance  of  ^ from  different  wavenumbers, 
.■->pectra  such  as  those  in  figure  24  were  used  only  for  this  purpose  --  to 
determine  wavenumbers  (or  frequencies)  where  the  noise  dominated  over 
the  signal,  and  to  determine  the  upper  limit  of  the  sums  of  the  Fournier 
coefficients  of  equation  B-I4b.  For  this  visual  display,  only  the  relative 
contributions  are  needed  from  the  various  frequencies,  and  the  values  of 
the  ordinate  in  figure  24  are  in  arbitrary  units. 

2 2 

-ipcctra  of  (Uj)  , (u^) 

The  energy  density  spectra 


B-21 


Using  the  expressions  B-l6  to  B-18 

cp.,,(k,)=  ojIcIllL 

(IJL  \l2jj  A/2  l250f>5,i/V 


(2IL  \(2jj  ) 

O 


//go  A//^tP  cn)c 

- U5o^5,i/*J  ja 


_ .IJO .2  CC3)C*Ci) 

{2SOf^S.i/^J  ,2 


\ •=  U 2 » 


3=^ 


.V-/ 


B-2. 
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The  log(pj^(k2)  log  Cp22(1^3)  versus  log  spectra  are  plotted 
in  figures  18,  19,  2Z  and  23.  The  plotted  values  have  been  band  averaged 
to  give  2 average  values  per  octave.  The  value  of  C(i)^*(i)  was  compu- 
ted for  each  transform  value,  and  averaged  over  the  frequency  band.  The 
first  four  Fourier  transform  values  were  plotted  individually,  and  were 
not  averaged  together.  Next  the  5th  and  6th  transform  values  were  ave- 
raged, as  were  the  7th  and  8th.  The  next  band  was  comprised  of  the  9th, 
10th,  11th,  12th  values,  and  the  next  band  included  the  13th  to  1 6th  trans- 
form values,  etc. 

The  corresponding  wavenumber  is  the  geometric  mean  of  the  en- 
tire interval,  which  produces  the  minimum  distortion  on  the  log  k plots. 
For  example,  the  arithmetic  mean  of  C(i )C*(i ) for  j=9  to  12  formed  one 

j ' 

band  average.  The  average  of  these  was  plotted  at  the  wavenumber 


k = 


(8.5 


X 12.5/^--^^ 


B-23 


To  determine  the  fit  of  the  universal  curve  to  the  ^^nd 

(p22^^3^  spectra,  the  point  jlog  V , log  is  marked  on  the  log  ^Piilk^) 
or  log  vs.  log  k^  graph,  and  a line  of  slope  +1  is  drawn  through 

the  point.  The  point  (0,0)  on  the  log  F2(k/kg)  vs.  log  k/k^  graph  will 
lie  on  this  line.  Its  position  is  determined  by  the  fit  of  log  F^  to  the  log 
^ I j or  log  <p22  curve;  its  absissa  is  the  point  logC  ''^^  ^ and  the  ordinate 

•/  V 

is  the  point  logt"'^"  on  the  log  <Pj  ^ vs.  log  kj  graph.  From  these  points, 
the  \’alue  of  e^can  be  calculated.  A full  discussion  of  this  procedure  is 


given  by  Stewart  and  Clrant  (1962). 
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Universal  Dissipation  Spectra  G2(k/kg) 
From  section  2.4 

or 


(k^'f cp„(ki)  - (^sf  Gz&^ks) 


where 


2.  25 
2.  23 


ITT  J 


^3)  v^vC^O 
andcPjjjik^  ) is  defined  in  equation  B-22. 


B-18 


The  refore 


r /fc/  1 f Cfj)c'‘(j) 

{jzK  'ksj  L ^=.O^S.V'*  / 

^ifto  a V CCj)C*(j) 

~ ^ 250/35, v2  J 


250pS,v  J 


j= 

2 


B-24 


and  these  values  have  been  plotted  in  figure  20.  The  dissipation  £ was 
determined  from  the  fit  of  the  log  j j and  log  <P  curves  to  the  universal 
curve. 


145 

APPENDIX  C 

Signal  Handling 

(a)  Linear  gain  preamplifier 

Electronic  circuits  prior  to  July  12,  1974  utilized  a high  impedance 
preamplifier  located  near  the  shear  probe,  and  a differentiating  amplifier 
in  the  main  body  of  the  Camel.  These  circuits  are  shown  in  figure  C-1.  j 

The  shear  probe  has  been  represented  as  an  a.  c.  voltage  generator  in 
series  with  a capacitor  Cp.  Calibrations  performed  on  an  early  calibrator 
allowed  the  probe  and  preamplifier  to  be  calibrated  together.  The  calibra- 
tor described  in  Appendix  A,  completed  at  the  beginning  of  the  ATLANTIS 
II  cruise,  required  a preamplifier  to  be  wired  permanently  into  it,  and  a 
new  but  identical  preamplifier  was  built  for  this  purpose.  The  calibration 
preamplifier  and  the  field  preamplifier  were  identical  in  design,  and  wort' 
similar  in  gain  to  within  about  4%. 

The  output  voltage  from  the  preamplifier  is 

Wc,  ~ ^ fS.Vu.  A-4 


A (’ronmpl  i f ier 
H D i ffprcnt iat inp  Amplifier 


Sc 

1 .5xl0"^'^f 

\ 

lO^^H 

‘■2 

1 

C 

X 

R., 

2 . 

Txin^A 

c.,, 

1.8x10"^  f 

s 

1 . 

6x1  Jl 

ICl 

ADSnp, 

R7 

1 

xlo4  JI, 

TC2 

Mil 

s 

1 

xin^U. 

ri-jMre  f.-l  The  calibration  and  field  preamplifier  located  near  the  probe 
and  llic  d i ffcrf!n1:  ini  iiifi  am|'lificr  located  in  the  main  body  of 
1 he  Camel  in  use  for  profile  18. 
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and  again  a similar  expression  holds  for  the  output  voltage  lToz  of  channel 
2.  The  amplifier  differentiates  the  input  signal  to  give 


-JcolTci  = - dUci  = -JL  pS.V  jU' 

^ dt 


C-4 


and  also  attenuates  the  high  frequencies.  The  voltage  from  the  amplifier 

1. 414 

was  stepped  down  by  2.  5 by  the  FM  transmission  system,  so  the  output 
voltage  to  the  Brush  recorder  and  the  digitizer  is 


riMi4] 

1 “RjCz  1 

f Rt  + Rs  1 

r 1 eS.Vciu..'] 

L 2 sJ 

Q + jtaR3C3X' 

L 2 d-t  i 

C-5 

d tti  _ 

•2  5 •]rci  + j^*^3C3)(i+ 

f 1 

[ ■'  1 

L Rr+RgJ 

If  Taylor's  hypothesis,  valid  in  this  case  is  applied  (equation  3.  3) 

LToi 

C-6 

The  gain  and  high  frequency  attenuation  were  determined  by  meas- 
uring input  and  output  voltages  at  various  frequencies.  A good  fit  is  pro- 
vided by  the  values: 


_ -I339  (^1  +.00l7j‘<^)('‘'’  °0'Sj^)  lToi 
<5^3  pS,V^ 

c^u.  _ -2  08  0 -oo  i7jto)(  I + 00(8  lToz 
for  channels  1 and  2 respectively. 


C-7a 


C-7b 


(b)  Differentiating  preamplifier 

The  pyroelectric  effect,  noted  in  section  3 generated  a low  fre- 
quency voltage  in  the  probe  sufficiently  large  to  saturate  the  preamplifier 
fluring  free-fall  profiles  at  the  equator.  To  reduce  the  low  frequency  vol- 
tages in  the  preamplifier,  a special  differentiating  preamplifier  was  con- 


structed  for  use  during  free-fall  profiles  called  a field  preamplifier,  and 
a band  pass  amplifier  with  uniform  gain  over  the  frequencies  of  interest 


was  used.  The  calibration  preamplifier  was  unchanged.  These  circuits 
are  shown  in  figure  C-2. 

The  output  of  the  field  preamplifier 


C-8 


(jpi  _ -jcoR.F Cp 

( I + jcoRiF  Cif)  . 

and  the  output  of  the  calibration  preamplifier, 

I r _ lT; 

- C.c 

may  be  combined  to  eliminate  the  capacitance  of  the  shear  probe  Cp  which 


is  an  unknown. 


J-p,  ^ jC^RirC,c 


+Rt 


LTc, 


_ RiF  Cic 

C I + JLu>R|FCir)  . 


R&+  Rt 


•’  aSyi  dui 

s ^ -1  dT 


The  response  of  the  amplifier  is 


C-9 


lT  = 


-'/c 


L 


ju>  RiCj. 


4 1 


r 


4 J 


r lTfi 


C-10 


and  ayain  the  FM  system  steps  down  the  voltage  by  1.414.  When  Taylor's 

2.  5 

hypothesis  is  invoked,  the  final  expression  for  the  shear  is 


du.  .,[2.5  -ir 
L\.4ihJl 


2.5 


- R 


’/R  = 


L 1?,FC,r  . 


T?f 


LTo  > 


+ c-11 


fo r channe 1 1 . 

The  three  circuits  (two  preamplifiers  and  the  amplifiers)  were 
tested  by  measuring  the  input  and  output  voltages  at  various  frequencies 
while  on  board  ship,  and  also  by  a w'hitc  noise  test  utilizing  the  PD  P- 1 2 
compviter  at  the  Institute  of  Oceanography.  Eiach  channel  of  each  individual 
circuit  was  tested  separately:  the  probe  was  simulated  by  an  a.  < . power 
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Thn  cal  ihral  inn  and  differentiating  fielil  prcampl  i f' iers 
located  near  the  shear  probe;  and  the  amp  1 i 1' ie'r  located  in  the  ■ 
body  of  the  Camel,  in  nse  for  pi'ofile  '.’n  nnd  I'ol  lowing  prid'iles 


